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Kearneysville, West Virginia 25430 





Dear Colleague: 


The Eastern Energy and Land Use Team (EELUT) is pleased to provide you 
this report on the results of a regional survey of the chemistry of 
headwater lakes and streams in the New England states. This study was 
conducted to assess the status of these waters with respect to acidifi- 
cation. The report is a part of the EELUT series on air pollution and 
acid rain. Other reports in this series are listed on the inside front 
cover. 


The extent and location of surface waters that are acidified or vulner- 
able to acidification was determined from a survey of 226 headwater 
lakes and low order streams in six New Snqland states. This data was 
compared to historical temporal trends, and were used to evaluate pro- 
posed models of acidification and to identify potential indicators of 
acidification. A companion study (FWS/0BS-80/40.17) provides similar 
information for the middle Atlantic states. 


Please teel free to send suggestions or comments to £ELUT so that we 
may continually strive to improve our future products. 


Sincerely, 
i+ 2 {7 ‘ . 
fi. Kes} Sehtek« ‘ 


x. Kent Schreiber 
Acting Team Leader, EELUT 
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EXECUTIVE SUMMARY 


A survey was conducted of 226 headwater lakes and low order streams 
in the six New England states. The waters selected were relatively 
undisturbed by direct human disturbance and were low in color. The sites 
were stratified by geographic area, and by bedrock geology and soil 
cation exchange capacity. 


Acidic (pH <5) surface waters occurred in every state. 
Approximately 8% of the waters surveyed were pH <5, and 29% were pH <6. 
The low pH waters tended to occur in clusters, although high pH waters 
frequently occurred in the same area. Waters with low alkalinity content 
were more common than low pH waters, and tended to occur in the same 
areas as the acidic waters. Approximately 24% of the waters had 
alkalinity cogcentrations of 20 eq 1 or less, 41% had concentrations 
of 100 weq 1 “ or less, and 53% were 200 ueq 1 “ or less. The 
proportion of acidic surface waters in New England is similar to that in 
other areas where bedrock is low in buffering capacity and precipitation 
is similarly acidic. Hydrogen ion content was not correlated with color, 
an index of organic acids. Therefore the acidic lakes are not primarily 
the result of natural organic acids. Highly colored, acidic lakes were 
excluded from the survey. 


Calcite saturation index (CSI) results were similar to alkalinity, 
and CSI and alkalinity were highly correlated. Approximately 59% of the 
waters surveyed had CSI values greater than 3, and were classed as 
susceptible to acidification. This is the same conclusion reached besed 
on alkalinity data. Specific conductance, an index of total ionic 
concentration, was correlated with alkalinity and could be used as 2 
simple index of vulnerability to acidification. However, the coefficient 
of correlation was low. 


Te lakes and streams surveyed contained appreciable sulfate 
concentrations (80-120 yeq 1 “) that could not be attributed to marine 
aerosols. The concentrations were similar to those found in other areas 
where precipitation is similarly acidic but were relatively uniform, and 
sulfate was not positively correlated with hydrogen ion as observed in 
other surveys. The lack of correlation may result because precipitation, 
and thus sulfate deposition, is chemically uniform over the region 
surveyed. 
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Aluminum and manganese concentrations of the lakes and streams 
surveyed were correlated with hydrogen ion contert, although the 
correlation was weak for manganese. The regression for aluminum was very 
similar to that obtained from other regional surveys and apparently 
represents the thermodynamic equilibrium between hydrogen ion and 
aluminum., Of the 226 waters surveyed, 13 had aluminum concentrations of 
200 ug 1 or more and pH of 5.5 or lower, conditions that may be toxic 
to sensitive fish species. 


Of the physical factors measured, the ones most highly related to 
buffering capacity and acidity were bedrock geology and, to a lesser 
extent, soil cation exchange capacity. A sensitive lake or stream, as 
determined by alkalinity, was nine times more likely to be found in an 
area where bedrock was low in bufferinro capacity. A sensitive lake or 
stream was three times more likely to be found in an area where soils 
were low in cation exchange capacity. Soil class was correlated with 
bedrock class, with sensitive soil classes laraely found in areas where 
bedrock was also sensitive. In headwater areas soils tend to be thin and 
therefore may be less important than bedrock in these areas in 
determining water chemistry. 


The only other physical factors of importance were lake area and 
stream order. Lakes of all sizes were low in alkalinity but oly lakes 
of 20 ha area or smaller were acidic. Higher order streams were higher 
in pH and alkalinity. This suggests that some factor related to 
watershed size may be an important factor in neutralizing acid 
deposition, and watershed size was correlated with lake pH for lakes in 
Maine. The waters surveyed were selected to be low in human disturbance. 
Disturbance was a!so correlated with bedrock geology, probably as a 
result of topography, confounding the effects of disturbance on water 
chemistry. 


Historical pH and alkalinity data indicate that waters located in 
areas where buffering capacity is low have been acidified. Of the 95 
lekes with usable historical pH data, 61 (64%) were lower in pH in this 
survey. Of 56 lakes for which historical alkalinity data were located, 
39 (70%) were lower in the present survey. The average lake for which 
historice! data were available increased in hydrogen ion content 
five-fold and decreased in alkalinity by 60%. 


The relationship between pH and calcium content of surface waters 
(Henriksen model) was tested as an indicator of acidification. This 
model predicted that 57% of the waters surveyed had been acidified, a 
value in ar agreement with that obtained from lakes for which 
historical water chemistry data were found. The model further predicted 
an average loss of about 50 »peq 1 © of alkalinity. Historical date, 
all from low alkalinity lakes, indicated a decline of about 109 
veq 1". Given the larce range of chemical conditions encountered, 
this was considered to be reasonably good agreement. The Kramer mode) 
gave good agreement with the data, 75% of the samples were within 0.5 pH 
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unit of the predicted pH based on carbonic acid weathering alone. The 
remaining 25% of the samples were more than 0.5 pH unit lower than 
predicted and are presumed to be acidified. The Thompson mode! also gave 
good agreement with the data and a similar prediction of acidification. 
The Kramer and Thompson models yield similar predictions of acidification 
that are about half that of the Henriksen model. 


We conclude that a substantial portion of the headweter lakes and 
low order streams in New England are vulnerable to acidification, snd 
that alkalinity is the best measure of vulnerability. Bedrock geclogy 
was the best physical factor that could be used to predict surface water 
alkalinity, and thus vulnerability to acidification, but a substantia! 
portion of the waters that were predicted to be vulnerable were not. 
Both historical water chemistry comparisons and acidification models 
indicate thet the vulnerable lakes have declined in alkalinity and 
increased in hydrogen ion content, presumably as the result of 
atmospheric deposition of acid. 
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INTRODUCTION 


The precipitation that falls on the northeastern United States is 
highly acidic compared to that in much of western North America. The 
annual weighted mean pH of this precipitation ranges from 4.2 to 4.4 
(National Atmospheric Deposition Program 1981). idic precipitation has 
been identified as the cause of acidification of surface waters, and 
consequently adverse effects on aquatic oraanisms, in Scandinavia (Almer 
et al. 1978; Muniz and Leivestad 1980), Ontario, Canada (Harvey 1980), 
and New York (Pfeiffer and Festa 1980). 


The effects of acidic precipitation on aquatic ecosystems depends on 
a variety of climatic, a th hy topographic, morphometric, biotic, and 
anthropogenic factors (Galloway and Cowling 1978). The relative 
importance of these factors is presently poorly understood, and the 
identification and quantification of aquatic resources at risk from 
acidic precipitation has been inadequate. In the United States the 
effects of acidic precipitetion on surface water chemistry and aquatic 
biota has been documented only for the Adirondack mountain region of New 
York. Several studies have predicted that substantial portions of the 
northeastern United States are vulnerable to acidification based on 
indirect evidence (Gallowey and Cowling 1978; Hendrey et al. 1980), but 
these predictions have not been widely tested. 


A number of models of acidification have been proposed (Henriksen 
1980; Thompson 1982). These models predict vulnerability of surface 
waters to acidification besed on chemical interrelationships, and 
estimate the amount of acidification that has already occurred. 
models have not been tested in the northeastern United States. Existing 
water chemistry data bases are not satisfactory for this purpose as most 
sites are located in productive lowland areas, consist of large lakes or 
high order streams, and have been affected by direct anthropogenic 
activitity (e.g., agriculture, waste disposal) 


This study was conducted to collect appropriate water chemistry data 
to assess the chemical status of headwater lakes and streams in the 
northeastern United States with respect to acidification. From these 
data the extent and location of surface waters that are acidified or 
vulnerable to acidification was determined. The data were compared to 
historical temporal trends, and were used to evaluate proposed models of 
acidification and to identify potentia! indicators of acidification. 
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METHODS 
Selection of Sampling Sites 





Surface waters suitable for use in this study were selected from al) 
available waters according to the following: 


© geographic coverage 

o bedrock buffering capacity 

0 soil cation exchange capacity 

0 location in watershed 

o direct human disturbance in watershed 
¢ availability of historical data 

oO access 


The waters sampled were primarily lekes (193 of 226; 85%). Streams were 
sampled only when a suitable lake could not be located in a cesired area, 
or was not accessible either physically or legally. 


Geographic distribution of lakes over the region was as uniform as 
possible. A few areas were inadequately represented. For exemple, it 
was very difficult to locate undisturbed lakes with public access in 
Connecticut. 


Sampling sites were selected to represent al] classes of bedrock 
buffering capacity and soil cation exchange capacity present in the 
region. Bedrock buffering capacity classes were those described in 
Hendrey et al. + , and soil cation exchange capacity groups were 
those of McFee (1980). The original maps depicting the distribution of 
these factors were obtained from the authors ‘S. Norton, Dept. of 
Geology, University of Maine, Orono, Maine 04469; W. McFee, Dept. of 
Agronomy, Purdue University, West Lafayette, Indiana 47907). These maps 
were generally of 1:509,000 scale, the same as the U.S. Geologica! Survey 
state base map. Lakes were selected to represent the various categories 
in rough proportion to the area of the categories in each state. Lakes 
may have been misclassified because of errors in precisely locating the 
lake on the state base map (many sma]! lakes are not shown on these maps 
and locations were transferred from 74 or 15 minute quadrangle maps), or 
because cf loss of resolution on the geology and soil maps as a result of 
scale or smoothing. 
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Headwater sites were generally low in direct human disturbance of 
watershed. Many otherwise suitable lakes were rejected because of 
disturbance in the watershed. In some cases, if no alternate lekes or 
streams better fit the criteria, disturbed waters were sampled. An 
arbitrary code was devised that represented differing levels of 
disturbance encountered, and waters were assigned a code based on a 
Teele hia of the lake shoreline and visible portions of the watershed 

Table 1). 


In each state contact was made with the appropriate state agency 
responsible for water chemistry data records (e.9., fish and game 
department, water cuality department), and any federal entity managing 
aquatic resources in the state (e.g., U.S. Forest Service). Agencies 
were asked to supply existing historical water chemistry data as wel! as 
te suggest potential samp'ing locations. Particular attention war yiven 
to selectinc sampling sites for which there were historical! data 
available; about half (95 of 226) of the final sites had such data. 


Finally, accessibility cf the sites was considered in making a final 
selection. Preferred sites were at least ore half mile from the nearest 
vehicle access point, but not more than two miles from that point. Most 
sites had a marked foottrail leading from the access point to the lake. 


Althcvoh size was not a selection criterion lakes fitting the above 
criteria were cenerally sma!!. Lakes lerger than 75 ha were not sampled. 
Data on watershed area were available for lakes located in the state of 
Maine (M, Hutchins, Environmental Studies Center, University of Maine, 
Orono, Maine). Wetershed area was included as a variable in this subset 
of the data. 


Sample Collection Procedur’ 





Each water body selected for survey was visited by a two person 
survey team, Most waters were reached on foot, but a few were reached 
by fixed-wing aircraft or portoon-equipped helicopter. A standardized 
sampling procedure was followed: 


o Survey watershed: one person walked as much of the shoreline 
as feasible, taking note of the major vegetation, exposed 
bedrock, inflow and outf'ow streams, aquatic veoetation, and 
evidence of human activity (roads, campsites, seasonal 
dwellings, permanent dwellings, logging, agriculture). 

Such observations were made on the trip into and out of 
the area as well. 

o Collect water samples: one person waded to the center of a 
stream or rowed to the deepest area of a lake in an 
inflatable reft. In streams, water samples were collected 
from mid-depth directly into the sample bottles. In lakes, 
one set of samples was collected at 0.1 m below the surface 
directly into the sample bottles and one set near the 








Table 1. Criteria Used to Estimate Relative Amount 
of Human Disturbance 





~ Disturbance 





No visible disturbance: access 
by foot trail or aircraft only. 


Slicht disturbance: access 

by 4-wheei drive vehicle: one 

or two seasonal dwellings in 

the watershed; evidence of 

logging activity but not presently 
active. 


Moderate disturbance: access by 
2-whee! drive vehicle; more than 
two seasonal dwellinas; evidence 
of active logging but not 
immediately adjacent to water. 


Severe disturbance: access by 
paved road; permanent homes around 
lake; hn Saad in watershed; 
active loaqging adjacent to water. 











bottom (without disturbing sediment) in a plastic Van Dorn-type 
water sampler. Water temperature was measured with 2 pocket 
thermometer at the time of collection. Water samples were 
placed in ecid-washed, distilled water rinsed linear 
polyethylene bottles. Each set of samples consisted of three 
bottles -- one 500 ml for pH, alkalinity, specific conductance, 
and color; one 250 ml for sulfate and chloride; and one 125 m! 
for cations. The 125 ml bottle cortained 6.25 ml of 4N 

nitric acid. Each bottle (except the cation sample bottle) 
was rinsed with sample water twice before filling. Samples 
were placed on ice within 1 hour after collection. 


The samples were generally collected during the open water period, 
but a few were collected through the ice. Previous experience indicated 
that, except for the spring run-off period, water chemistry in these 
waters is relatively stable throughout the year. In any event, 
geographic variability far exceeds temporal variation in individual 
waters. The lakes sampled were not euthropic or dystrophic. Inasmuch as 
the surface and deep samples were chemically similar, there was no 
evidence that algal production or hypolimnion decomposition significantly 
affected water chemistry in these lakes. 


AVALYTICAL METHODS 
Field Procedures 





Analyses of pil, alkalinity, specific corductance, and color were 
performed at field locations. Within 8 hours after sampling, and as soon 
as possible, the 500 m! bottle was removed from ice and warmed to room 
temperature. Two 100 wm) aliquots were removed for determination of p# 
and alkalinity. The pH was measured with a portable meter (Fisher mode! 
107 or Cole Parmer DiaiSense) equipped with a plastic-body, gel filled 
combination electrode. The meter was standardized with pH 7 and 4 
buffers, and electrode response wat verified by measuring the pH of 
dilute sulfuric acid solutions of theoretical pH 4. If measured values 
deviated from expected values by more than 0.1 pH units the electrode 
was discarded. The electrode was rinsed are with dis -illed 
water, blotted dry, and soaked in the sample for 15 minutes or longer -- 
until three successive readings at 1 minute intervals were identical -- 
and pH was recorded. 


Alkalinity was determined by titrating each of the 100 ml sample 
aliquots with 0.0200 N sulfuric acid to pH» 4. Acid was added in 
0.10 ml portions using a micro syringe until pH 5 was reached, then in 
0.05 ml portions to pH <4. The pH was recorded after equilibration 
following each addition of acid. Alkalinity was calculated by two 
methods. Inflection point alkalinity was determined by the method of 
Gran (Stumm and Morgan 1970), ard fixed endpoint (pH 4.5) alkalinity was 
determined as described in American Public Pealth Association et al. 
(1975). Inflection point results were used for all analyses and 
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comparisons except for those using historical data, where fixed endpoint 
data were used. 


Two 50 ml aliquots of sample were measured and used for 
determination of specific conductance and color. Specific conductance 
was measured with a calibrated meter (Markson Scientific Company model 
10), and apparent color was determined by comparison cf unfiltered 
samples with platinum cobalt standard solution (LaMotte Chemical Company, 
Chestertown, Maryland). 


Laboratory Procedures 





The remaining water samples were kept on ice, returned to the 
laboratory, and kept refrigerated until analyzed. The unacidified 
250 ml water sample wes used for measurement of sulfate by the 
microthorin method (Fritz and Yamamura 1955), and chloride by the 
ferricyanide method (American Public Health Association et al. 1975). 
All analyses were performed in replicate. The acidified 125 ml sample 
was used for measurement of cations. Sodium and potassium were measured 
with air-acetylene flame atomic absorption spectrophotometry (AAS; Perkin 
Elmer model 703), calcium and magnesium by nitrous oxide-acetylene flame 
AAS, and manganese by graphite furnace AAS. Aluminum was measured 
initially by the colorimetric method of Dougan and Wilson (1974), and 
later by graphite furnace AAS. The two methods gave comparable results. 
As samples were unfiltered, the results obtained represent total 
concentrations of the various ions. 


Calcite saturation index (CSI) was calculated from the equation 
given by Kramer (1976). For samples where alkalinity was zero or 
negative, a modified equation (Kramer 1981) that approximates CSI was 
used, 


A subset of 31 samples ranging in color from 0 to 260 units were 
analyzed for total organic carbon (TOC) by G. Glass, U.S. Environmental 
Protection Agency, Duluth, MN. A linear regression analysis was 
performed to determine the relationship between color and TOC. 


Quality Assurance 





The functioning of analytical instruments was assured by frequent 
analysis of known standards. Response of pH, specific conductance, 
spectrophotometers and atomic absorption spectrophotometers was checked 
daily. Most analyses were performed in replicate to enable calculation 
of variance. U.S. Environmental Protection Agency Water Pollution 
Control Samples for Trace Metals, and Minerals (EPA Environmental 
Monitoring and Support Laboratory. Cincinnati, Ohio), were analyzed 
following the same procedures used for water samples. An interlaboratory 
exchange of samples was conducted with the Norwegian Institute for Water 
Research. 
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As a further check on analytical eccuracy and data coding errors, we 
calculated tcte! fonic balance and theoretical specific conductance. Al! 
jon concentrations were converted to units of .eq 1 “, anions and 
cations were summed, and the sum of anions was compared to the sum of 
cations for each semple. Theoretical specific conductance was 
calculated by multiplying fon concentrations by equivalent conductance 
values for 25°C obtained from the Handbook of Chemistry and Physics 
(Weast 1978). Theoretical conductance was compared to measured 
conductance for each sample. 


Data Analysis 


Results obtained from the above analyses were entered into a SAS 
(Statistical Analysis System, Raleigh, NC) data set operating or an IBM 
model 4341 computer located at the University of Maine, Orono, Maine. 
SAS statistical analysis and graphic routines were used for data 
analysis. Punch cards produced from this data set were sent to 
Paul Rago, U.S. Fish ard Wildlife Service, Ann Arbor, MI, where they 
were entered into a MIDAS (Michiaan Interactive Data Analysis System) 
data set using the University of Michigan computer system. This system 
was employed for analysis of variance and multiveriate data analysis. 





Historica! Data 





Historical pH and alkalinity data were obtained through 
correspondence with water chemistry and fisheries offices in the six 
states. All such data were for lakes. Data were screened for obvious 
errors and only data for which the method of collection was documented 
and reasonable were incorporated in this study. 








RESULTS AND DISCUSSION 


Quality Assurance 





The laboratory analyses of EPA lity Assurance Samples for 
Minerals (Table 2) and Trace Metals (Table 3) cave acceptable accuracy. 
The coefficient of variation was always less than 5% for minerals and 
less than 10% for trace elements. The element determined with poorest 
precision was aluminum; however, the recoveries we obtained were 
comparable to those reported by EPA, and our standard deviations were 
smaller. An interlaboratory comparison was conducted with the Norwegian 
Institute for Water Research for calcium, magnesium, sodium, and 
potassium determinetions on 2] water samples collected in Norway, and 
sulfate and chloride determinations on 18 water samples collected in 
Maine. The following linear regression equations were obtained: 

Calcium (ME) = 1.02 Calcium (NIVA) -0.C3 

Meonesium (ME) = 0.97 Magnesium (NIVA) + 0.003 

Sodium (ME) = 1.06 Sodium (NIVA) + 0.01 

Potassium (ME) = 0.97 Potassium (NIVA) = 0.03 

Sulfate (ME) = 0.94 Sulfate (NIVA) + 0.33 

Chioride (ME) = 1.09 Chloride (NIVA) -0.19 
In all cases intercepts were not different from zero and slopes were rot 
different from one (t test, p <0.05). 


Recovery of known additions of trace metals and minere!s was always 
within 5% of that expected. The pH of dilute sulfuric acid solutions was 
always within 0.1 units of theoretical. In those instances where results 
deviated from theoretical by more thar this amount, the electrode was 
replaced and the meter recelibrated. Acceptable results were then 
obtained. 


A further check on the accuracy of water chemistry determinations 
was made by comparison of the sum of anions and the sum of cations 
(Figure 1), and of theoretical and measured conductance (Figure 2). The 
sums of ions were nearly equal, with only two obvious outliers. The 
line of best fit had an intercept of 14 veq 1 = and a slope of 1.00, 
not significantly different from 0 end 1.0 res ryt (p 0.0001). 

The linear correlation was highly significant (r° = 0.97, p 0.0001). 
The line of best fit for theoretical and measured conductances had an 
intercept of -1.22 wS cm “, and a slope of 1.12, both not 
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Table 2. Results of Analysis of EPA Water Pollution Quality Control Samples for Minerals. 
Mean Values for PH were Computed from Hydrogen lon Concentrations 
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EPA sample True _ Laboratory results 
number Factor value vT 3.0. tv. or TT. 
(N=3) 
? pH 7,4 7.53 ° 2.3 7.51-7.57 

(units) 
Calcium 6.7 6.8 0.2 7.9 6.5-7.1 
(mq/1) 

Magnes ium ?.4 2,4 0.05 2.1 2.32.5 
{mq/1) 

Potassium 1.7 1.6 0.06 3.8 1.5-1.7 
(mg/1) 
Sodium 7.0 }.? 0.06 0.8 7.1-7.3 
(mq/1) 

4 pH 8.6 8.56 - 4.6 £.51-8.67 

(wnits) 
Calcium 32.0 4.8 0.? 0.6 4. 3-35.39 
(mg/1) 

Magnes ium 7.1 7.1 0.1 1.4 6.9-7.3 
(mq/1) 

Potassium 1.2 7.5 0.06 0.8 7,4-7.6 
{mq/1) 

Sodium 40.0 40.1 0.1 0.3 79, 9-40, 3 
(ma/1) 


bis! . - burl 











Al 


Al 


Table 3. 

Trwe 

Value ' 
(Ne) 

350 474.3 

55 67.3 

50 68.3 

il 13,7 

70 776.7 

350 w7.5 
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Figure 1s) Sum of total enions versus sum of tote! cations 
in veq 1 ~ for 196 sites, surface samples only. Least 


squares regression eouation: 
Cations (r° = 0.97), 
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Fioure 2. Theoretical versus measured specific conductance 

for 195 sites, surface samples only. Leest squares 
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significantly different from 0 and 1.0, respectively (p@.0001). The 
linear regression was highly significant (r~ = 0.97; p = <0.0001). 
These two comparisons indicate that there are no major analytical or 
coding errors in the water chemistry data set. 


Site Distribution 





we sampled 226 lakes and streams in the six New England states 
(Figure 3). Many other sites were visited and rejected as unsuitable for 
inclusion in this study, because of high color, disturbance in the 
watershed, and other reasons. The distribution of sites among states 
{Table 4) was roughly proportional to the state area. Lakes comprised 
85% of the sites and streams 15% (Table 4). Streams were sampled only 
when a suitable lake could not be located in the desired area. The lakes 
sampled were either seepage lakes with no visible inlet (N=122) or 
drainage jlakes with permanent inlets and outlets (N=71). Stream order 
ranced from 1 to 4 (N: 1215, 2*5, 3*2, 4=*11). 


The sample sites were distributed among bedrock geology classes in 
proportion to the area of each state that was underlain by each bedrock 
class (Table 5). Sites were selected in areas within each soil class 
present in each state (Tabie 5), but the proportional distribution was 
not determined. Most bedrock types were represented in rough proportion 
to the proportion of land area in each state, especially classes 1 and 
2. All soil classes present were represented, but the proportional 
distribution by land area is not known. 


CHEMICAL FACTORS 
pH and Alkalinity 





Surface waters of pH <5 were found in all six states (Figure 4). 
Nineteen sites (8%) were in this category. Sixty-five sites (29%) had 
pH <6. The low pH waters tend to occur in clusters, i.e. southeastern 
Maine, northern and southern Vermont, etc. A similar distribution was 
found for iow alkalinity waters (Figure 5). There were 53 sites (24%) 
with alkalinity <20 .eq 1°, 93 (41%) with alkalinity <100 veq 1 “, 
and 120 (53%) with alkalinity <200 veq 1 ~. 


A number of other ional surveys of surface water chemistry have 
been conducted. Results from the distribution of lakes by pH range for 
areas where precipitation is acidic (pH<4.6 weighted annual average) and 
non-acidic (pH>4.6) are summarized in Table 6. There is a pronounced 
difference in the proportion of acidic lakes (pH<5) between areas that 
do and do not receive acidic precipitation. Within the acidic 
precipitation group, however, there is considerable v‘riation in the 
proportion of lakes in each pH range. This is not surprising 
considering the number of factors that may affect lake pH (e.g., organic 
acids, bedrock geoloay, soil chemistry, wastershed size, land use, 
etc.), and the variability in the procedures followed in the different 
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Figure 3. Map of New England showing the location of the 
276 sample sites and type of water body sampled. 
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Table 4. Number of Lakes and Streams Sampled in Each State 








State Number of lakes Number of streams 
Connecticut 16 7 

Maine 62 0 
Massachusetts 34 0 

New Hampshire 38 1] 

Rhode !sland 8 0 
Vermont _ 35 _ 15 

Total Number 193 33 

Grand Total 226 
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Table 5. Distritutton of Semple Sites Among Bedrock Geology and Soi! Class. Bedrock Geology 
Classes are those Described by Hendrey et 21, (1990), and Soll Classes are those 
of Mefee (1980). A Dash Indicates that thet Category was not Present 











funter of Percent of state ares Number (and *) of sample sites Sumber (and ¢) of sample 
sites in bedrock cless in bedrock class sites in soil class 

1 ? T re a | — SS SS 
Commer t tcwt 73 * ws 5 2 l 71 1 0 1% 4 6 
(4) (9?) (4) ($7) (17) (26) 
Maine 6. 1? 44 20 x we 1m 6 0 7 B) 2? 
(79) (61) {1) (it) ($3) (%) 
Massachusetts i4 Me 31 27 10 16 6 6 ? ? 18 + 
(%) = (19) {?) (6) {?1) ($3) (76) 
few Hampshire 49 w 6K ’ ] 1] 4H 0 0 17 4 
(22) (7) (35) (65) 
Phede 1s land i 6) % ? 0 j ] 0 . 3 ? 3 
(fa) (1?) (#8) (24) (38) 
Vervunt 4) } 6) w 16 § y i 4 ? 4h 
(1) (&) (16) (@) (4) (9) 
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Figure 4, Map of New Encland showing the distribution of 
surface water samp’es by pH. Date are for surface samples 
only and are means of two replicates. 
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Figure 5. Map of New England showing the cistribution of 
sample sites by alkalinity. [fata are for surface samples 
only and are mears of two replicates. 
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Table 6. Regional Water Chemistry Survey Results for 
Surface Water PH Distribution 








Location Number Percent in pH range Reference 





| 





Areas Where Precipitation Averages <pH 4.6 


New England 226 g 21 71 This Study 

West Sweden 314 36 21 43 Aimer et al. 1974 

West Sweden 15 27 47 27 Dickson 1975 

South Sweden 51 2 20 78 Malmer 1975 

South Norway 155 18 38 44 Wright et al. 1977 
South Norway 719 64 33 3 Wright and Snekvik 1978 
Denmark 14 29 57 14 Rebsdorf 1980 

Scotland 72 26 36 38 Wright et at. 1980 

Nova Scotia 2] 52 24 24 Watt et al. 1979 

Quebec 25 12 40 48 Jones et al. 1980 
Central Ontario 26 8 58 34 Scheider et al. 1979 

La Cloche Mountains 152 28 34 38 Beamish and Harvey 1972 
Sudbury 150 13 15 72 Conroy et al. 1976 
Adirondack Mountains 849 25 30 45 Pfeiffer and Festa 1980 


Areas Where Precipitation Averages »pH 4.6 


North Norway 77 0 13 87 Wright and Gjessing 1976 
Northwest Wisconsin 265 0 6 94 Lillie and Mason 1980 
North Minnesota 85 0 0 100 Glass and Loucks 1980 
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surveys. For example, Wright et al. (1977) followed a randomization 
procedure to select an unbiased set of lakes in south Norway, whereas 
Wright and Snekvik (1978) working in the same area surveyed only lakes 
for which fish population data were available, although they avoided 
very large or highly disturbed lakes. New England falls within the low 
end of the range of the proportion of acidic lakes, being similar to 
south Sweden and central Ontario. Although we selected small, 
headwater, undisturbed waters for our survey, ye attempted to obtain a 
representative sample from a large (168,000 km”) area with wide 
chemical diversity. Many other surveys were restricted to areas where 
buffering capacity was known to be low. 


Regional surveys that included alkalinity measurements are less 
common (Table 7). Again, there is a pronounced difference between 
regions where precipitation is acidic and regions where it is non-acidic 
with respect to the proportion of low alkelinity waters found. Waters 
with total alkalinity <20 neq 1 ~ or <100 weq 1 ~ are much more 
common in areas that receive acidic precipitation. New England is 
intermediate in alkalinity in the group of regions that receive acidic 
precipitation. As the regions are all known to be low in buffering 
capacity, the increased proportion of very low alkalinity waters in acid 
rain districts may be the result of titration of alkalinity by deposition 
of excess hydrogen ion. 


Calcite Saturation Index 





CSI has been proposed as an index of vulnerability to acidification 
that is superior to a simple measure of buffering capacity (Kramer 
1976). Generally, larger values of CSI reflect lower buffering 
capacity. We grouped our samples into categories of CSI <1, 1-3, and »3 
(Figure 6.) These categories were modified from Kramer (1976) and Glass 
and Loucks (1980). CSI <1 represents waters that are saturated or 
nearly saturated with respect to calcite and would not be susceptible to 
reduction in pH from acidic precipitation. Values of 1-3 represent 
waters that are potentially susceptible to acidification, and values of 
>3 are waters that are susceptible to acidification. Of the waters 
surveyed, 15% are not susceptible, 26% are potentially susceptible, and 
59% are susceptible. The distribution of susceptible waters is very 
similar to that obtained from pH and alkalinity results. This is not 
surprising because CSI is calculated from both pH and alkalinity, as 
well as calcium. CSI was highly correlated with ajkalinity (Figure 7). 
The regress ten equation was log alkaiinity (weg 1 “) = 3,42-0.39 CSI 

.929, p <0.0001). The majority (59%) of the waters surveyed in 
New England were classed as sensitive to acidification, CSI >3. This 
proportion is very similar to that classed as sensitive on the basis of 
alkalinity alone (53% were alkalinity <200 ueq 1 “). 

Regional surveys that included measurement of Calcite Saturation 
Index are rare. Kramer (1981) calculated CSI with water chemistry data 
collected at road crossings over streams by e variety of workers in 


20 





Table 7. Regional Water Chemistry Survey Results for Surface Water 


Alkalinity Distribution. 


Alkalinity values are in veq | 





Location Number 


< 


Percent in alkalinity ranae Reference 


> 





Areas Where Precipitation Averages <pH 4.6 


New England 226 23 
South Norway 62 3 
Denmark 14 86 
Nova Scotia é} 71 


Central Ontario 26 

Ontario 600 2 
La Cloche Mountains 4 100 
Adirondack Mourtains 692 4) 
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14 


12 
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15 


18 


47 
pA 


16 


This Study 

Wright et al. 1977 
Rebsdorf 1980 

Watt et al. 1979 
Scheider et al. 1979 
i)cmerman and Harvey 1979 
Beamish 1976 

Pfeiffer and Festa 1980 


Sreas Where Precipitation Averages »pH 4.6 


Northwest Wisconsin 265 - 
North Minnesota 85 0 


15° 


17 
26 


68 
52 


Lillie and Mason 1980 
Glass end Loucks 1980 
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b £1 - 200 veg l 
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Figure €, Map of New England showing the distrButioneof 
sample sites by Calcite Saturation Index (CSI). Data are 


for surface samples only, and are mears of two replicates. 
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Figure 7. Pelationship of alkalinity to Calcite Saturation 
Tnoex for 195 New England lakes and streams, Date are means 
cf two replicates and means of surface and deep samples 
where taken. Least squares negression equation: 

log Alk = 3.42 - 0.39 CSI (r° = 0.929). 
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eastern Canada. The data were presented as a contour map, and showed 
that portions of central and western Ontario had CSI values of 3 or 
more, and that most of Quebec and the Maritime Provinces had water with 
CSI greater than 3. Glass and Loucks (1980) calculated CSI for 85 lakes 
in northern Minnesota, an area that does not receive acidic 
precipitation. They found 11% had CSI <1, 53° had CSI of 1-3, and 36° 
were CSI>3. The distribution of New England headwater lakes enc 

streams with respect to CSI indicates that many more are sensitive to 
acidification than in Minnesota (59° vs. 36%), and that the distribution 
may be similar to that in Ontario. However, CSI does not seem to provide 
any superior predictive ability to pH and alkalinity measurements, and 
represents an additional mathematica!) manipulation. Zimmerman and Harvey 
(1979) reached a similar conclusion for surface waters in Ontario. 


lonic Composition 





Specific conductance is an index of total ionic content of water, 
and thus has potentia! for use as a simple index of buffering capacity. 
However other factors, such as distance from seacoast and differences in 
ion activities, also affect specific conductance, limiting its 
usefulness for this purpose. As a result, the relationship between 
specific conductance and alkalinity (Figure 8), although strtistically 
significant, is not especially strong (r° = 0.71). Some waters that 
are high in specific conductance are low in alka'inity. In some cases 
this represents increased hydrogen ion content, which "3s very high 
equivalent conductance (Weast 1978). In other cases it represents 
marine aerosol input, or unusual minerology such as gypsum. 


Chloride was measured for use in correction of other ions for marine 
aeroso! contribution, and is a possible indicator of human disturbance. 
The relation of chloride to distance from the sea coast (shortest 
straight line) is given in Figure 9. The influence of marine aerosols 
close to the coast is evident, with high chloride concentrations close to 
the coast falling off very sharply at, 75 to 100 km to background 
concentrations of less than 50 veq 1 “. This is almost identical! to 
the situation in south Norway (Wright and Henriksen 1978), and Scotland 
(Wright and Henriksen 1980). There are four waters with anomalously high 
chloride concentrations. These four waters also were ranked high in 
human disturbance. Mowever, there were many other sites that were also 
ranked high in human disturbance that could not be identified on the 
basis of high chloride concentration. As chloride concentration was 
almost entirely related to marine aerosols, it was used to correct other 
ions for marine aerosols contribution. Chloride was not used as an index 
of human disturbance. 


lonic composition of the surface waters, corrected for marine 
aerosols, varies markedly with pH (Figure 10). Total fonic content is 
highest at high pH and declines with successively lower pH. Of the 
anions, bicarbonate decreases markedly below p¥ € and is absent below 
pH 5, as expected. However, although sulfate composes a larger 
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Figure 8. Relationship between alkalinity and specific 
conductance for '°5 New Enqland lakes and streams. fata 
are means of two replicates and means of surface and deep 
csemples where taken. Least (es 0 we ccuation: 
Pik = -~185 + &.0CE Sp. Cond, *« 0.711) 
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Figure 9. Relationship of chloride concentration to 
distance from seawater for 184 New England lakes. 
Chioride data are means of two replicates and means 
of surface and deep samples where taken. Disterce 
from seacoast is the shortest straioht line to open 
ocear. 
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Figure 10. Distribution of major ions in New Encland 
lakes and streams clessed according to ph renge. 
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proportion of total anions as pH declines, the concentration cf sulfate 
does not change appreciably. Of the cations, hyd fon increases as 
pk declines, as expected, calcium and magnesium decline mertedly, end 
sodium and potassium decline moderately. The present composition of 
cations and anions is compared to world averages in Table 8. Hich pH 
New Enoland lakes have anion distribution similar to the world aversces, 
but @ somewhat 'eroer proportion of ca'cium and magnesium. 


The mean tote) tonic content of the surface waters surveyed wes 
about 1,000 »eq 1°° at pH»6, 300 veq 1°° at pH 5-6, and 200 veq 1 
at pH<S (Figure 10). Henriksen (1980) summarized for composition data 
for 13 surveys of arees unaffected by acidic precipitation. Mean total 
jonic content ranged from 200 to 2,300 eg 1°, but im 11 of the 13 
surveys the mean was less than 800 veq 1 “. There was no apparent 
cifference in tonic strength of acidic New England lakes and non-acicic 
softwater lakes elsewhere in the world. 


Sulfate 


Virtually a7] of the sulfate in the New England lakes and streams 
surveyed is from non-marine sources. One lake, Holts Pond, 
Massachusetts, hec an extremely high sulfate copcentration \mean of 
surface and deep samples, replicated = 44 mg 1° 915 weg } ). The 
sulfate was nearly belarced by calcium (mean 16.3 mg 1 *, 815 
ee 1°), indicating that this lake wes most likely irfluenced by a 
local deposit of gypsum, and it was net included in sulfate analysis. 

For the remaining waters, the SO,: C! ratio (Figure 11) shows that most 
samples exceeded the ratio of these fons in seawater (0.10), especially 
those locatec more than 75 km from the seacoast where marine eerosols are 
relatively unimportant. This indicates a non-marine source of sulfate. 

As further evidence of this, tote! sulfate concentration is independent 
of chloride concentration (Figure 12). The S0,: Cl ratios are higher 
than those reported for south Norway. Wright dna Henriksen (1978) report 
S0,: Cl ratios 3-20 times higher than in seawater. In the lakes we 
sufveyed the ratios were as much as 70 times higher than seawater, and 
40°. of the lake’ ' ratios exceeding 70 times thet in seawater. 


The conce: +e*’ ns of total and non-marine sulfate found in 
New England lekes «a: « comparable to those found in other rections where 
precipitation is acidic, and are much higher than in areas where 
precipitation is not acidic (Table 9). However, contrary to the 
situation in south Norway, non-marine sulfate cortent was not positively 
correlated with hydrogen fon content in waters with low (<100,eq 1°") 
calcium content (Ficure 13). Brown and Sadler (1981), using the date of 
Wrioht and Snekvik {197e}, found a highly significant correlation in 471 
lakes in south Norway (r° = 0.354, p 0.0001), even though they did 
not stratify the lakes by calcium concentration. Brown et 1. (1980), 
using the same data, also found significant positive correlations 
between non-marine sulfate and calcium, enc calcium plus magnesium, in 
accition to hydrogen jon. These correlations were low in New England, 








Table 8. Percent Relative Distribution of Major Ions in New England 
Headwater Lakes and Streams as Compared with World 
Averages, Corrected for Marine Aerosols. Ion 

Concentrations are veq 1” 











Tons world averages” New England 

Clark Livingstone Rodhe pHs pH 5-6 pH ¢ 
Cations 1924 1963 1949 — — ; 
Ca + Mg 88.2 67.9 88 .6 79.5 19.9 91.0 
Na + K 11.8 12.1 11.4 20.5 20.1 9.0 
Anions 
HCO. 83.0 82.0 83.0 0 15.4 74,° 
50, 17.0 18.0 17.0 100 84 .6 25.4 





"After Henriksen (1980) 
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Figure 1l. Sulfate: chloride ratio versus distance 
from sercoest for 184 New England lakes. Sulfate 

and chloride are means of two replicates and means 

of surface and deep samples where taken. [Cistance 
from seacoast is shortest straight line to open ocear. 
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Figure 12. Relationship of sulfate concentration to 
chloride concentration. Sulfate is total concentration, 
not corrected for marine aerosols. Date are for 184 New 
England lakes, are means of two replicates, and are 
means of surface and deep samples where taken. 
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Table 9. Amount of Sulfate (veq 1°*) Found in Surface 
waters in Recional Water Chemistry Surveys 





Region ‘ Tota!) sulfate WNon-marine sul fete Reference 
7 3.0. Tt 3.0. 








Areas where Precipitation Averages pH 4.6 


New England 22 138 $1 27 ae This Study 

South Yorway »€ 99 33 92 - Sjessing et al. i976 
west Sweden é 200 70 18c ° Grahn et al. 1974 

Scot land 7? 178 9) 48 76 wright and Giessing 1976 
Nowa Scotia 16 1$2 2€ 11? 25 watt et al. 1979 

Central Ontario 15 | SA 97 - - Dillon et a]. 1978 
Sudbury, Ontarto ‘ £04 292 800 290 Scheider et al. 1975 

ia Clocne Mayntainre 8 2a? 8: as 4) Beamish 1976 


Areas Where Precipitation Averages py 4.6 


Central Norway 6? 6§ 34 63 - Glessing et al. 1976 
western Ontartc ] 63 . cP . Armstrong and 


Schindler 197] 














120 





. 
100 
H 860 
I 
6 
N 
60 
U 0 0 
; 
/ 4o o 
L 
0 .) 
20 '.) ol lCUP o 
O 
oO . 
o «USP o 0 
0 oe) . 
0 0 o0 @ Gann® G9? cote SP w 
wore qeevwuwewe ere qoeewewewe Goverowuuuuurereeegeeeeee ee egyeyeteeee weerre 


' ' | 
0 20 40 60 80 100 120 140 =: 160 
NON-MARINE SULFATE (UEQ/L) 


Figure 13. Relationship of hydrogen ion te non-mar ine 
sulfate concentration for 72 Kew England waters with 
calcium <100 veq 1 “. Data are means of two replicates 
and means of surface and deep samples where saken. 

The Hers 9 regression was not significant (r° = 0.047, 
p = 0.12). 
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r= = 0.096 and 0.070, respectively. In our samples, high calcium 
Zeters had significant bicarbonate ion whereas low calcium waters were 
lacking in bicarbonate and sulfate was the major anion. The lack of 
correlation between surface water pH and non-marine sulfate content may 
result from the fact that precipitation throughout the region is 
relatively uniform in chemistry (National Atmospheric Deposition Program 
1981) unlike Scandinavia, except for the relationship between marine 
aerosols and distance from the seacoast. Waters receive relatively 
uniform amounts of non-marine sulfate deposition, but only some lakes 
become acidic as the result of low buffering capacity. 


Another possible explanation is that there is a regiona! 
“background” level of sulfate, in addition to marine aerosols, that masks 
the sulfate associated with acid. Background sulfate could originate 
from terrestria! dusts or anthropogenic sulfate that is not acidic. 

Wright (Richard Wright, Norwegian Institute of Water Research, personal 
communication) has proposed that a proportion of non-marine sulfate 
equivalent to the sum of camera sodium plus potassium be subtracted 
as a correction for background, sulfate. Background sulfate calculated by 
this method averages 50 veq/! © for New England. Background sulfate 
content of a lake should be independent of ionic strength if it is truly 
of regional character. However because of the principal of ionic balance 
there will probably be some relationship between background sulfate and 
total ionic strength. Background sulfate calculated from our data is 
correlated with total gations (Figure 14), but the slope is very close to 
zero. Using 50 veq 1 “ as the best estimate of regional background 
sulfate, net non-marine sulfate was calculated and compared to hydrogen 
ion (Figure 15). The relationship is no better than that obtained with 
non-marine sulfate. 


Calcium and Magnesium 





9 Calcium and magnesium were highly correlated with alkalinity 
(r“ = 0.951, p <0.0001; Figure 16). The regression soeatie -{¥ 
Alkalinity (eq 1*) = -55 + 0.95 (non-marine Ca + ly 
This compares favorably with rersauagion Alkalinity ney 1” ' = «14 + 
0.93 (non-marine Ca + Mg, veg 1 = 0.98 obtained by Henriksen 
(1980) based on 13 sets of data from areas that are relatively 
unaffected by acidic precipitation. The slopes are very similar but the 
intercept is lower for our data. This suggests that acidic 
precipitation has not affected the proportional relationship of calcium 
plus magnesium to alkalinity in New England, ee alkalinity may be 
reduced generally. In contrast to this, Almer et (1978) reported 
that the relationship between these factors was nearly 1:1 in an area of 
Sweden that received a smal] sulfur load, but much less than 1:1 in an 
area that received a higher sulfur loading. Sulfur loadi 
relatively uniform over New England (National Rtmsoheric fepesition 
Program 1981). 


Acidic lakes tend to be lower in calcium than comparable non-acidic 
lakes. Almer et al. (1974) found that calcium plus magnesium content 
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Figure 14, Relation of calculated background sul fete 
concentretion to the sum of catiors for 168 lakes. Lakes 

with non-marine sul fete or calcium plus magnesium concentration 
exceeding 400 vec 1° were eliminated. The regression equation 
is BSU = 0.26 cations - 3.49, r° = 0.52, p = 0.0001, computed 
from means of two rep!icates and means of surface and deep 
samples where taken. 
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Figure 15. Relation of hydrogen jon to net non-marine sulfate 
concentration... The regression equation is # ion = 3.24 + 0.02 
net sulfate, r~ = 0.006, not significant. 
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Figure 16. Relationship of alkelinity to non-marine 
calcium plus magnesium for 226 surface waters in 
New England. Data are means of two replicates and 
means of surface and deep samples where taken. 
Least squares regression equation: 

Alk = ~-53 + 0.93 (CatMg) (r° = 0.948). 
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increased from 290 veq 1"! in Swedish lakes of pH 4.9 or less to 820 
yeq 1 ° in lakes with pH 7.0 or more. Gjessing et a]. (1976) reported 
that lakes ip south Norway with median pH 4.76 had mean calcium content 
of 57 yeq 1" whereas lakes in central Norway with median pH 6.40 had 
mean calcium content of 139 veq 1 “. The explanation for this is that 
areas that receive acidic precipitation contain both well buffered and 
poorly buffered waters, depending on the chemical content of bedrock and 
soil. Well buffered waters contain proportionally larger amounts of 
calcium and bicarbonate, which is of geological origin, than poorly 
buffered waters. Only the poorly buffered waters will be acidified, and 
therefore acidic waters will be lower in calcium than non-acidic waters 
in the same area. Mohn et a]. (1980) concluded that headwater lakes in 
south Norway that were low Tn calcium, and therefore in buffering 
capacity, had water concentrations of hydrogen and sulfate ion that were 
significantly related to precipitation noe poe al Precipitation 
chemistry in the New England region is relatively uniform, and this 
relationship cannot be tested. 


Aluminum and Manganese 





Aluminum concentration was related to (Figure 17), with high 
aluminum concentrations generally found in, Ow pH waters, as expected. 
The regression equation is log.) Al (ug 1°) g 3.93 - 0.36 pH and 
the correlation coefficient is Rignificant (r° = 0.37, p 0.001). -] 
Eighteen of 226 samples (8%) had alumi gum concentrations of 200 yg 1 
or higher, the highest being 510 ug 1 ©. These waters were 
distributed in the same areas as low pH waters (Figure,18). Forty-seven 
samples (21%) had aluminum concentrations of 100 yg 1 “ or more. 
There, were 13 sample sites that had an aluminum concentration of 200 
vg 1 ° or more and a pH of 5.5 or less. These conditions may be toxic 
to some species of fish (Schofield and Trojnar 1980). 


The relationship between aluminum and pH found in New England is 
very similar to that found in lakes in Sweden (Dickson 1975, 1980), 
Norwey (Wright and wpe Raed Scotland (Wright and Henriksen 1980), 
and New York (Schofield 1982), as shown in Table 10. Schofield and 
Trojnar (1980) made multiple measurements of pH and aluminum 
concentration in one stream in New York and found a similar relationship 
to that found in the regional surveys cited above (Table 10). This 
suggests that the relationship is an expression of the thermodynamic 
equilibrium of aluminum and hydrogen ion, enabling the prediction of 
aluminum concentration from pH in any region. 


Manganese is only slightly related to pH (Figure 19). Manganese 
generally was found in lower concentrations than aluminum, and there is 
more variability in the manganese: pH relationship (r° = 0.015; 

p »0.05). The reason for this variability is unknown. 


Organic Acids 





Color was measured as an index of organic acids. Color was highly 
correlated with total organic carbon (TOC) determinations performed on a 
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Figure 17. Relationship of total aluminum concentration 
to pH. The pH data are means of two replicates, and both 
pH and aluminum are means of surface deep samples where 
taken. Least squares regression equation: 

log Al = 3.93 - 0.36 pH (r* = 0.37). 
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Map of New England showing locatiog of sample 
sites grouped by aluminum concentration (ug 1°). Data are 
for surface samples orly. 
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Table 10. Values for the Intercept (A) and Slope (8) for the Linear 








Regression of Total Aluminum and pH. Log,, {ug 1 “) = A +8 (pH), 
from Regional Water Chemisty Surveys 
Location N Intercept Slope Reference 
New England 226 3.93 -0.% This Study 
Sweden 37 3.85 -0.35 Dickson 1975 
Sweden 73 3,94 -0.38 Dickson 1980 
Norway 47 3.85 -0.41 Wright and Gjessing 1976 
New York 214 4.14 -0.33 Schofield 1982 
New York 68° 3.97 -0.87 Schofield and Trojnar 1980 





"Acid Brook, measured 68 times from March to August. 
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Figure 19. Relationship of manganese to pF for 195 
New England surface waters. The pl data are means of 
two replicates, ar. both pH and mangenese are means 
of surface and dees samples where teken. Least squaress 
recression equation: e 

log Mn = 1.83 - 0.05 pH (r* = 0.015). 
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set of 31 samples —- im color from 0 to 260 units. The regression 
ry was TOC (mg “Rs -32 + 0.0613 color (color unit), and 

This equation R. used to calculate TOC values for al! 
This relationship agrees well with the approximation TOC = 
color : 10 used by the Norwegian Institute for Water Research 
(R. Wright, personal communication). 


A plot of color or log, TOC versus p¥ (Figures 20,21) shows no 
relationship between these tors. Altheowoh we attempted to avoid 
highly colored waters in the selection of sampling sites, a few were 
included. On some occasions water turned out to be more colored than we 
realized when the sample was collected. In other cases a colored water 
was sampled if there were no obvious bog characteristics and we needed a 

le in that area. The lack of relationship between pH and color ‘or 
TOC} ts that the acidic lakes are not the result of organic acids. 
Similar ty there was little relationship between alkalinity and corel 
(Figure 22). Only 16 samples had alkalinity of more than 500 veg | 
and these al! had color less than 4 units. The remaining 210 waters 
showed no relationship between alkalinity and color. There is no 
evidence of increased buffering capacity (alkalinity) contributed by 
organ acids in highly colored waters, as was reported by Johannessen 

) for humic lakes in Norway. 


Driscol! et al. (1980) reported that aluminum formed strong 
complexes with organic ligands. 't is therefore possible that changes in 
total organic carbon, which are independent of pH, could account for some 
of the variation in a) umi num concentration in the waters surveyed. A 
comparison of aluminum and total orgaric carbon concentrations (Figure 
23) indicates that this relationship is not significant (r® « 0, 7", 


Precipitation Chemistry 





To successfully correlate surface water and precipitation chemistry, 
a large — of highly uniform geology and climete with strong 
gradients in precipitation chemistry would be required. Such a situation 
exists in Scandinavia (Wright and Gjessing 1976), but not in New England. 
As a result, attempts to relate surface water chemistry to precipitation 
chemistry were unsuccessful. Chemistry of precipitation in New England, 
as reported by the National Atmospheric Deposition Program (1981), is 
relatively uniform. The small differences in annual deposition of 
hydrogen ion over the region are far outweighed by differences in surface 
water chemistry resulting from watershed interations, bedrock type, sot! 


type, etc. 
PHYSICAL FACTORS THAT AFFECT WATEP CHEMISTRY 


Elevation 





Surface water pH was only weakly related to elevation (Figure 24). 
Some very acidic lakes are located at elevations below 100 m, and lakes 
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Figure 20. Relationship of pH to apparent color. Data are 
means of two replicates and means of surface and deep samples 
where taken. 
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Figure 21. Relationship of total organic carbon (TOC) to pH. 
TOC was calculated from color by the equation: 
TOC = 1.32 + 0.0613 color. 
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Figure 22. Relationship of alkalinity to apparent color. 
Data are means of two replicates and means of surface and 


deep samples if taken. 
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Figure 23. Comparison of aluminum concentration to total 
organic carbon concentration in 196 New England lakes and 
streams. 
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Figure 24. Relationship of elevation of sample site to pH. 
Elevation was determined from U.S. Geological Survey topographic 
maps. The pH data are means of two replicates and means of surface 
and deep samples if taken. 
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of pH <5 are distributed throughout the range of 10 to 1,000 m 
elevation. However, waters with pH of 6.5 or more were seldom found 
above 600 m. A similar situation exists with alkalinity (Figure 25). 
Surface waters with negligible alkalinity were found at al! elevations 
from 10 to 1,000 m, but all but two, lakes lociuted at 600 m elevation or 
higher had alkalinity of 100 »eq 1 “ or less. 


The influence of elevation on surface water chemistry is most likely 
a complex interrelationship of bedrock geology, soi! type and depth, 
depth of glacial till, and watershed size. Higher elevation areas would 
be expected to have bedrock that is resistant to chemical weathering, and 
thin soils and till. Watersheds may be smaller at higher elevation 
because of the steepness of the terrain. The presence of low pH and low 
alkalinity waters at low elevation indicates that marine influence does 
not necessarily result in increased buffering capacity in New England, in 
contrast to the situation in Norway. Some very acidic lakes are located 
near the seacoast in Maine, Massachusetts, and Rhode Island. In Norway 
acidic lakes were rare at elevations less than 200 m, because of the 
presence of readily weathered sediments deposited during the late 
Pleistocene marine submergence (Wright and Snekvik 1978). Although the 
proportion of acidic lakes and streams is greater at higher elevations in 
- England, there are as many acidic lakes at lower as at higher 
elevations. 


Size 


In general, acidic lakes were also smal! lakes. Al! but two lakes 
that had a pH of 5.5 or less were also less than 20 ha surface area 
(Figure 26). Obviously, not al! small lakes were acidic, but lakes 
larger than 20 ha rarely were. But the larger lakes were low in 
alkalinity (Figure 27). All but lakes with area greater than 20 ha 
also had alkalinity of 400 veq 1°* or less. Other regional surveys of 
lake chemistry seem not to have considered lake area as a factor 
governing chemistry. In general, larger lakes would be expected to have 
larger, more complex watersheds with a greater variety of bedrock types 
and thicker soils. Therefore we expected larger lakes to be better 
buffered and higher in pH (lower in hydrogen ion). The largest lakes we 
surveyed were relatively low in alkalinity, and therefore buffering 
capacity, yet tended to be of higher pH than smaller lakes of comparable 
alkalinity. 


There is apparently some factor related to size that caused the 
larger lakes to have lower concentrations of hydrogen jon than smaller 
lakes of the same alkalinity. One possible explanation is watershed 
size. Larger lakes tend to have larger watersheds, which may be capable 
of neutralizing more hydrogen ion even though they are low in bufering 
capacity. This could be the result of lower-gradient slopes and thicker 
soils and till, which would allow more contact time. We tested this 
hypothesis using lakes in Maine, for which drainage area data were 
available. We compared 10919 drainage area to pH, drainage area to 
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Figure 25. Relationship of sample site elevation to alkalinity. 
Elevation was obtained from U.S. Geological Survey topographic 
maps. Alkalinity data are means of two replicates and means of 
surface and deep samples where taken. 
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Figure 26. Relationship of lake area to pH. The pH data are 
means of two replicates end means of surface and deep samples 
where taken. 
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Figure 27. Relationship of lake area to alkalinity. Alkalinity 
data are means of two replicate and means of surface and deep 
samples where taken. 








alkalinity, and the ratio lake area/drainage area to pH and alkalinity by 
use of linear regression. The only ey eowy significant regression 
was for 10949 drainage area and pH (r° = 0.14, p <0.003; Figure 28). 


This suppo the hypothesis that increased oes of the drainage basin 
can neutralize increased amounts of hy¢rogen ion independent o 
alkalinity. 
Drainage Type 

Lakes were classed as drainage lakes (with outlet) and s lakes 


(no outlet). The two types of hydrology spperentty do not significantly 
affect lake water chemistry (Table 11). Orainage lakes hed a nearly 
significant tendency to be more highly colored than seepage lakes. In 
contrast to this, lake hydrology was an important factor in controlling 
+ i? capacity in northern Wisconsin lakes. (Eilers 

et al., in preparation). 


Lakes and streams were chemically similar, differing only in pH, 
sodium, and potassium (Table 12). The streams were higher in pH, which 
may reflect the acid neutralizing effect of laroer watershed area or the 
fact that lakes can receive an appreciable amount of precipitation 
directly on the lake surface, without watershed contact. 


Stream Order 


Higher order streams tended to be 4 in pH and alkalinity, 
specific conductance, Calcite Saturation Index, calcium, magnesium, 
sodium, potassium, and chloride (Table 13). There were no siginifcant 
differences for aluminum, manganese, sulfate, or color. This pattern is 
as expected, because higher order streams have larger watersheds and are 
more likely to include e variety of bedrock geology types. 


Bedrock and Soil Class 


Bedrock logy type had a pronounced effect on chemical factors 
related to acidity and ee capacity, but little effect on other 
factors (Table 14). Mean pH, alkalinity, specific conductance, calcium, 
and magnesium increased consistently along the spectrum from bedrock A 
l ‘ low buffering capacity) to type 4 (very high buffering capacity 

cite saturation index consistently decreased. 








Waters located in areas where bedrock was class 1 or 2 varied widely 
in pH (Figure 29) covering nearly the entire range measured. Waters 
underlain by bedrock class 3 or 4, however, were largely restricted to 
higher pH. Low alkalinity waters were much more abundant in areas with 
bedrock class 1 and 2, and high alkalinity waters were more abundant in 
class 3 and 4 (Figure 30). However, low alkalinity waters were found in 
all bedrock classes. 


The bedrock geology classes devised by Hendrey et a1. (1980) were 
intended to reflect buffering capacity inherent in the bedrock. this. 
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Figure 28. Relationship of 109) drainage area to labe pH. 
Lake pM values are means of two licates and means of 
surface and deep samples where tater. Least squares 
regression equation: 

log drainage area ~ 6.26 + 0.31 (pH). 
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Table 11. Mean Values and Standard Deviations of Chemical Factors 
for Lakes Classed by Hydrology Type. P is Probability that 
Means are Drawn From the Same Populations, One Way 
Analysis of Variance Test 


Factor Hydrol t p 
Draina Seepage 
i Eo. $.0. 








PH, units 6.36 0.86 6.27 0.91 0.66 
Alkalinity, veq 17! 231 411220483 0.88 
Specific conductance, uS ca”! 54 Ss 49 46 0.51 
Calcite saturation index, units 3.67 1.72 3.52 1.67 0.62 
Calcium, mg 17! 5.4 ee ee ee 
Magnesium, mg 17! 0.9 0.7 Ll 19 0.39 
Sodium, mg 17 2.3 2.6 2.5 3.0 0.61 
Potassium, mg 1”! 0.6 0.5 0.6 O58 0.95 
Aluminum, vg 17! 85 111 a9 113 0.83 
Manganese, vg 1”! 40 46 38 30 0.80 
Chloride, mg 17! 3.6 5.2 3.5 4.6 0.95 
Sulfate, mg 17 4 5.7 06.4220 
Color 30 28 23 19 = s«0.06 
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Table 12. Chemica! Factor Means and Standard Deviations for Lakes and 
Streams. P “s Probability that Means are Drawn From the Same 
Populations, One Way Analysis of Variance Test 














Lake Stream P 
¥ “$3.0. ¥ 3.0. 

PH, units 6.29 0.89 6.78 0.53 0.01" 
Alkalinity, veq 17 281 534 217 130 0.55 
Specific conductance, uS cm $4 53 57 28 0.81 
Calcium, mg 17 5.4 8.4 4.0 2.0 0.40 
Magnesium, mg 17! 1.1 1.7 1.2 0.7 0.73 
Aluminum, vg 17! B4 101 1 82 0.55 
Manganese, .g ~ 39 35 45 4&5 0.44 
Sodium, mg 17! 2.2 2.6 3.5 2.5 0.02" 
Potassium, mg 1” 0.6 0.5 1.2 0.5 0.01" 
Chloride, mg 17! 3.2 4.6 6.2 5.1 0.06 
Sulfate, mg 1” 6.5 2.5 6.9 2.3 0.46 
Color, units 14 80 12 66 0,39 





"significant fora = 0,05 





Table 13. Relationship of Stream Order to Water Chemistry Factors. 
Values Reported are Arithmetic Means of a)! Available Samples. 
Means with the Same Superscript Letter are not significantly 
Different (One Way Analysis of Variance Test 
with Scheffe's Comparison, p <0.05) 











Factor Stream order 

5 2 3 T 
PH, units 6.3" 6.2" 722° 7,0°% 
Alkalinity, weq 17 174° gst = ggz®eSod ogg A 
Specific conductance, wS oa”? 43° 35° 5q*s 716 
Calcite saturation index, units 3.5° 3.8° 2,5°°* 2.6?" 
Calcium, mg 17 en ae on 4.52 
Magnesium, mg 1”! 10° 0.5% 17° 13°" 
Sodium, mg 17 2.27 2,48 get 4.9° 
Potassium, mg 17! 0.9% 0.78 —s_ gto? 1 59 
Aluminum, vg 1” 127° 110° 14° a9* 
Manganese, vg 1” sa* 53° 13" 39° 
Chioride, mg 17! 25% gat g gael g.1° 
Sulfate, m 17! 6.0 —s «5, 5" 54° 7.6 
Color, units 35° 42° 10° 21° 
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Table 14. Mean Values of Chemical Factors for Surface Waters Classed 
by Bedrock Geology. Bedrock Classes: 1 = Very Low Buffering 
Capacity, 2 = Moderately Low Buffering Capacity, 3 = 
Moderately High Buffering Capacity, 4 = Very High 
Buffering Capacity (Hendrey et a1. 1980). Means 
with the Same Superscript Letter are not 
Sianificantly Different (One-Way 
Analysis of Variance with 
Scheffe's Comparison 











p <0.05) 
Factor Bedrock class 

I z 3 
PH, units 6.02" 6.25° 6.66" 7.21° 
Alkalinity, veq 17 g* 172° 575° ~—s-1., 025° 
Specific conductance, vS cm”? 3° 43° ga? 149° 
Calcium, mg 17? 2.1° 3.9° 11.5° 18.7° 
Magnesium, mg 17! 0.7? 0.8° 1.5° 4.8° 
Aluminum, ug 172 90° 80° 110° 114° 
Manganese, ug 17 41° 34° 38° 102° 
Sodium, mg 17* 2.8° 2.0° 31° 2.5° 
Potassium, me 1") 0.5° 0.6° 0.8° 0.8° 
Chloride, mg 17* 4.2° 3.0° 3,2° 2.8° 
sulfete, mg 17 5.5? 6.52 7.6? 15.9° 
Color, units 26° 26° 25° 28° 
Calcite saturation index, units 4.08" 3.87° 2.03° 1.68° 
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Figure 30. Frequency distribution of surface waters by bedrock 
geology in alkalinity intervals. Each alkalinity value given is 
the midpoint of the interval, + 50 veq 1 ~. 
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geological factor does affect surface water chemistry for those chemical 
factors that are related to acidity and buffering capacity but not other 
factors, as intended. The relationships are consistent but lack 
precision because of high variability within a bedrock class. Waters 
located in bedrock classes 1 and 2 were not significantly different from 
each other for any chemical factor measured, and classes 3 and 4 were 
different for only a few factors. Based on the above analysis bedrock 
types 1 and 2 were pooled and types 3 and 4 were pooled. The ability of 
these groupings to predict surface water sensitivity to acidification 
was tested by preparing a contingency table for bedrock geology and 
alkalinity. Bedrock typés 1 and 2 were termed sensitive and types, 3 and 
4 were termed non-sensitive. Waters with alkalinity of 290 yeq | * or 
less were termed sensitive and those exceeding 200 veq | “ were termed 
non-sensitive. The hypothesis that bedrock class has no effect on 
alkalinity was tested by use of a contingency table (Table 15). From 
this table we calculated ,°, the log-likelihood ratjo (G-statistic), 
and the odds ratio (Bliss 1967, Fleiss 1981). The x and G statistics 
were highly significant © See The odds ratio indicated that a 
sensitive lake or stream is nine times more likely to be found in bedrock 
class 1 or 2 as in class 3 or 4, 


There were only three of the five possible types of soil groups 
(McFee 1980) in New England: SS,, SS., and NS. The highly sensitive 
groups (S1, $2) were not present. Sof] type was much less important than 
bedrock type in affecting surface water chemistry (Table 16). Classes 
SS1 and SS2 were not significantly different for any factor measured. 

SS1 and SS2 were both significantly different from NS for alkalinity, 
calcium, sodium, and chloride. Class SS2 alone was different from NS for 
pH and sulfate. 


Soil classes in our sample were distributed nearly normally with 
be to pH (Figure 31), with intermediate pHs being most common in all 
soil classes and waters of all pHs being found in nearly every soil type. 
The relationship of soil classes to alkalinity (Figure 32) is highly 
skewed with an abundance low alkalinity waters. But the proportion of 
soi) classes within each alkalinity interval is nearly equal, and waters 
of all alkalinities are found in each soil class. 


The predictive ability of soil class for alkalinity, calculated in 
the same manner as for bedrock geology, is shown in Table 17. Soil 
classes SS1 and SS2 were combined as sensitive for this comparison. The 

and G statistics were highly significant ( xy p@.001, G p ©.005). 
The odds ratio indicates that sensitive lakes and streams are three times 
more likely to be found in areas with sensitive soi: class than in areas 
with non-sensitive soils. Thus soil class is not as good a predictor of 
alkalinity as is bedrock class. Kaplan et al. (1981), working with 
different bedrock and soil classes, concluded that there was no direct 
influence of bedrock on surface water chemistry but there was for soi! 
type, and bedrock indirectly influenced water chemistry through soils. 
However, the data set used in Kaplan's study consisted largely of 
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Table 15. Contingency Table Tests of the Predictive Ability 
of Bedrock Geology for Surface Water Alkalinity. 
Table Entries are Numbers of Lakes Sampled. 





Aikalinity class 











sensitive, Non-sens itive 
Bedrock class <200 veq 1 200 veq | Total 
Sensitive 
(Class 1 & 2) 135 29 164 
Non-sensitive 
(Class 3 & 4) 4a 18 27 
Tota! 144 47 191 





Test statistics 


a 29.98 (p <.001) 


" 


“4 


G 25.75 (p <.001) 


* 9.31 (p <.\\l) 


S.E. = 4.26 











Table 16. Mean Values for Chemical Factors Classed by Soil Cation 
Exchange Capacity. Soil Class: SS, = Slightly Sensitive Soils 
Dominate Area, SS, = Slightly Sensitive Soils are Significant 
but Less than 56% of the Area, NS = Mostly Non-Sensitive 
Soils (McFee 1980). Means with the Same Superscript 
Letter are not Significantly Different (One-Way 
Analysis of Variance with Scheffe's Comparison 











p <0.05) 
Factor Soil type 

3S, S, Ts 
PH, units 6.22%» 6.06" 6.46° 
Alkalinity, veq 17! 1S7e 83° 331° 
Specific conductance, vS cm”! 49° 4a? 55° 
Calcium, mg 1" 3.2° 2.5° 7,.0° 
Magnesium, mg 17! 1.1° 0.8° 11° 
Aluminum, vg 17! 83° 77° 95° 
Manganese, vg 17 35° 40° 40? 
Sodium, mg 17! 3.3° 3.1° 1,7° 
Potassium, mg 17 0.7 0.6" 0.5° 
Chloride, mg 17! 5.2° 5.6° 1,7° 
Sulfate, mg 17 6.12" 5.8" 7,8 


Color, units 27° 26° 25° 
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Figure 32, Frequency distribution of waters by soil class in 
alkalinity intervals. Each alkalinity value given is the 
midpoint of the interval, + 50 veq il. Class 3 = SS), 

4-« SSo. 5 = WS. 
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Table 17. Contingency Table Tests of the Predictive Ability of 
Soil Cation Exchange Capacity for Surface Wate. Alkalinity. 
Table Entries are Numbers of Lakes Sampled. 





Alkalinity class 
Sensitive “sensitive 











Soil class «200 weg 1~! 200 veq 1 Total 
Sensitive 
(SS, ‘ SS. ) Fo 14 102 
Non-sens itive 
(NS) 58 33 $1 
Tota! 146 47 193 





Test statistics 


© 13.26 (p -.001) 
G * 8.84 (p «.01) 

3.58 (p «.05) 
S.£. © 1,29 
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lowland sites, where soils would be thicker and more likely to influence 
water chemistry. In headwater areas soils are likely to be thin and may 
be absent in large areas, resulting in a greater influunce of bedrock on 
water chemistry. There my also be an interaction between bedrock class 
and soil class, with sensitive soils being associated with sensitive 
bedrock. Also, Kaplan et al. (1981) classified bedrock by formtion 
(igneous, metamorphic, etc.) rather than by chemical content. Thus, 
warble, for example, would be classed as metamorphic and relatively low 
in buffering capacity. 


Human Disturbance 





Human disturbence had an important effect on many water chemistry 
factors (Table 18). Generally, increasing human disturbance results in 
increased pH, alkalinity, calcium, magnesium, Calcite Saturation Index, 
and specific conductance. There was ro effect on aluminum, manganese, or 
color, and minimal effect on sodium and potassium. Chloride, often used 
as an index of human disturbance, di¢ not increase consistently with 
disturbance code but mean chloride values were higher in codes 3 
and 4 then in codes 1 and 2. 


Human activities such as logging, road and building construction, 
and agriculture disturb natural vegetation and expose soi! to increased 
leaching of chemical constitwents. Runoff may also be increased. 
Disposal! of human wastes further adds ions to the oe waters. 

This may explain the increese in chemical factors, especially those 
associated with the buffering complex, with increasing human disturbance. 
Alternatively, this relationsnip may result from an interaction between 
humen disturbance and bedrock geology. The waters surveyed in bedrock 
classes 1 and 2 tended to be low in disturbance whereas those in classes 
3 and 4 tended to be high (Table 19). This may reflect the tendency for 
human activity to be concentrated in areas where soi! fertility is high 
and topography is gentle, which are generally essociated with bedrock 
types high in carbonate minerals. 


The role of disturbance in the relationship among chemical factors 
was evaluated by performing stzet! tical analyses on the entire data set 
and on waters with disturbanc: < de i. or ? only. As there were no major 
discrepancies that could be 1’ *ributet to disturbance alone, we used al! 
sample sites in the analyses +» ported here. 


HISTORICAL COMPARISONS 
ph 


Of the 95 locations with wsable historical pH data, 34 (36%) either 
were the same or higher in the recent analysis and 61 (64%) were lower 
(Figure 33). If there were only random variation, it would be expected 
that the same proportion of locations would increase as decrease. The 
recent pH values are y ~{ Smo lower than historical values (paired ¢ 
test, t = 4.17, p = 








Table 18. Mean Values for Chemical Factors Classed by Disturbance Code. 
Disturbance Codes: 1 = No Disturbance, 2 = Slight Disturbance 
3 = Moderate Disturbance, 4 = Severe Disturbance. See Text 
for Description of Disturbance Assessment and Assignment 
of Classes. Means with the Same Superscript Letters 
are not Significantly Different (One-Way Analysis 
of Variance with Scheffe's Comparison, p <0.05) 

















Factor Disturbance class 

I 2 3 Z 
PH, units 6.10° 6.20° 6.88° 711° 
Alkalinity, yea 17! 138° 179° 412° 1,043" 
Specific conductance, »S$ cm” 33° 49° 81° 1319 
Calcium, mg 17 3.2° 4.5° 7,7> 14.6° 
Magnesium, mq 17 0.6° 0.9%» 1.5 5.0° 
Aluminum, ug 17! 82° 90° 79° 164° 
Manganese, ua 17) 37° 41° 46° 52° 
Sodium, mg 17! 1,3° 2.6° 45° 2,929 s¢ 
Potassium, mg 17! 0.4? 0.59 1.0° 1.0°°¢ 
Chloride, ma 17 1.9? 3.6° 6.0° 5.02% 
Sulfate, mg 17! 57° 7.3% 6.92 gg 4 
Color, units 23° 26° 31° 16° 
Caicite saturation index, units 4.0° 3,7° 2,7° 2.1° 
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Table 1%. 





Number of Surface Waters Sampled in Each Bedrock 


Class and Disturbance Code. See Text for Explanation 


ef Bedrock and Disturbance Code 














Disturbance Bedrock class 
Code I 2 3 
l 32 71 10 
2 47 98 16 
3 12 27 12 
4 2 0 4 
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Figure 33. Current versus historical pH for 95 New England 
lakes. Solid line indicates equivalent pH. 
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If the comparison is restricted to those waters that are sensitive 
to acidification, the proportion that have declined in pH is increased. 
For waters that presently have alkalinity of 100 neq 1° or less, 13 
of 53 locations (25%) either did not change or increased and 40 of 53 
(75%) decreased in pH (Figure 34). This differerce is also significant 
(paired t test, t = 4.38, p =<0.0001). A similar comparison for waters 
located in bedrock classes 1 or ? only (the lowest in buffering capacity) 
gave comparable results (Figure 35). 


The mean pH of stations for which historical data were located was 
5.37 (calculated from hydrogen ion) versus 6.07 historically. Or, stated 
another way, ean hydrogen ion increased from 0.8 ueq 1 * historically 
to 4.3 yveq 1 ~ recently, a five-fold increase. The historical data 
were determined as early as 1937 and as recently as 1978. Data more 
recent than 1978 were located but not used. There was no obvious trend 
of change in pH with elapsed time between measurements. However, the 
historical data set is rather small and includes lakes with a wide range 
in buffering capacity. Lakes with high buffering capacity would not be 
expected to change ry there ba in pH. If _gnly low buffering capacity 
lakes are considered (alkalinity <100 ueq 1 °) the data set becomes too 
small to be meaningful. 


The historical pH data comparisons, althouch not numerous, 
nevertheless indicate a substantial increase in hydrogen ion content of 
the lakes that are low in buffering capacity. The historical pH data 
were largely determined colorimetrically. Only one historical 
measurement was located that was determined electrometrica'ly. 
Colorimetric pH indicators may give erroneous results, especially in 
water with low buffering capacity (Boyd 1976, 1977, 1980; Pfeiffer and 
Festa 1980; Wright 1977). However, Schofield (1982) showed that the 
error was not large in Adirondack Mountain lakes, being 0.2 unit or 
less. We have done extensive testing of the colorimetric indicator 
methods commonly used in New Enaland and concluded that the limit of 
accuracy of the method was +0.2 unit, and the error was unbiased (Haines 
et al., in press). Norton et al. (1981) compared the two methods for a 
set of New England lakes and concluded that the two methods were 
equivalent within 0.1 unit. Almer et al. (1974) measured px in a number 
of lakes using both methods and found that the results were comparable, 
seldom differing by more than 0.2 unit. Lewis (1982) compared 
electrometric pH measurements with those made with the same type of 
indicators used 40 years previously and found “good agreement”. 


Many other comparisons of historica! and recent pH data have been 
made, A number of the most significant ones were reviewed by Haines 
(1981). Most comparisons that were made in areas where precipitation is 
acidic (pH <4.6 annual weighted mean) and where acid neutralizing 
capacity of bedrock is low have demonstrated similar declines in pH. 
Many studies in Scandinavia, reviewed by Wright and Gjessing (1976), 
reported that lakes declined from pH 6-7 during the period 1940-1950 to 
pH 5-6 in the early 1970's. Wright (1977) compared historical and 


7) 

















ee a ee 








 BEAAAAALASA BAAAABAADAS DAA AAASAAS LAB BBABE SE OG 


@ ~~ wo w => 
Voc cwze i+ a 


Figure 34. Current versus historical pH for 53 New England 
lakes with alkalinity of 100 uec 1 or less. Solid line 
indicates equivalent pH. 
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Figure 35. Current versus historical pH for 85 New England 
lakes located in bedrock class 1 or 2. Solid line indicates 


equivalent pH. 
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recent pH data for 128 lakes in southern Norway and found that 637 had 
become more acidic. Pecently, Morling (1981) reported annual pH 
measurements from 1960-1966 to 1979 for 10 Swedish lakes, with similar 
results. 


McColl (1981) found that the hydrogen ion activity of two reservoirs 
in the, Sierra Nevada mountains of California increased about 0.06 
yeq 1 during the period 1954 .o 1979. The mean monthly variation in 
hydrogen jon was correlated with mean monthly variation in air pollution 
in the San Francisco Bay Area. Lewis (1982) resurveyed 104 lakes in 
Colorado that had been surveyed previously. He found that pH declined 
fro a mean value of 7.09 during 1938-1960 to 6.87 recent: VMendrey et 
al. (1980) resurveyed 42 New Hampshire streams for which historical pH 
data were available and found that pH had declined in 37 (887). Ina 
previous survey of Maine lakes, Davis et a1. (1978) reported that the pH 
declined from 6.18 in 1937-1943 to 6.05 fm 1969-1974. This smal! decline 
most likely results from the fact that the lakes surveyed were large, 
lowland lakes that were not highly vulnerable to acidification. Pecently 
Norton et a). (1981) reported historical versus recent pH comparisons for 
69 lakes in Maine, New Hampshire, and Vermont and found that 51 (74%) had 
decreased pH. This finding agrees well with ovr results, and there is 
some overlap in the data sets. 


Alkslinity 





Historical alkalinity data were less common than pH. Only 56 lakes 
had such data. The historical alkalinity data were al! obtained by 
fived endpoint titrations to pH 4.5 (methy! orange endpoint). This 
method overestimates true alkalinity because the actual endpoint is not 
fixed but increased with oe alkalinity. We measured both fixed 
endpoint and inflection point alkalinity for our samples. A comparison 
of these two values (Figure 36) shows that they are highly linearly 
correlated (Inflection Point = -32 + 1.00 fixed endpoint; r° = 0.999). 
Thersfore fixed encpoint alkalinity data can be simply corrected for the 
overestimetion by subtracting 32 veq 1 ~, which was done for al! 
historical data. Conversely, we could have compared our fixed endpoint 
albalinity velues to unadjusted historical values. The results are the 
same in either case, but presentation of fixed endpoint alkalinity values 
indicates erroneously high alkalinity levels. 


A comparison of adjusted historical to recent alkalinity results 
(Figure 37) shows that 17 of 56 results (30%) increased or did not change 
and 39, (707) decreased. Historical alkaligity data averaced 166 
veq 1° and recent data averaged 68 »eq 1°, a 60% decrease. This 
decrease was significant (paired t test, t = 4.03, p = 0.002). 

Generally, waters with the highest historical alkalinity exhibited the 
largest declines. This may be because these waters had alkalinity that 
could be lost whereas waters that were low in alkalinity historically had 
little, to lose. The highest historical alkalinity value located was 600 
veq 1°, which is not very high on a world-wide basis. As for pH, there 
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FIXED ENDPOINT ALKALINITY 


Figure 36. ‘Irflection point (Gran plot) versus fixed endpoint 
(pH 4.5) alkelinity for 184 sample sites. Least squares 
regression equetion: IP = -32 + 1.00 FP (r° = 0.99%), 
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Fioure 37. Current versus historical alkalinity for 5é 
New England lakes. ‘Solid line indicates eouivalent alkalinity. 


























was no discernable relationship between time interval between 
measurements and magnitude of the change in alkalinity. The change in 
alkalinity acrees well with the change in pH and indicates comparable 
levels of acidification. 


Because historical alkalinity data are scarce, other comparisons of 
historical and recent alkalinity deta are rare. Morling (1981) reported 
annual alkalinity measurements for four Swedish lakes from 1967-1968 to 
1979. All declined by 10 to 20 eq 1°. Lewis (1982) reported that 
alkalinity averaged 440, eq 1 ~ in 64 Colorado lakes during 1938-1960, 
and averaged 36C yeq 1 © in 1979. McCol® (1981) analyzed weekly 
alkalinity data over the period 1944-1979 for one reservoir in the 
Sierra Nevada mountains of California. Alkalinity decreased from about 
350 veq 1 ~ to about 230 veq 1 ~. During the period 1954 to 1979, 
for which hydrogen ion measurements were also available, alkalinity 
declined about 90 ueq 1 ~. 


Multivariate Analyses 





The large number of factors measured in this study makes the 
identification of key variables difficult. A number of multivariate 
analytical methods were explored in an attempt to reduce the complexity 
of the data set and identify the factors that contribute most to the 
variance. 


Principal Components 





Principal Components Analysis (PCA) is a procedure that reduces a 
large set of intercorrelated variables to @ smaller set of hypothetical 
variables, each of which is independent of the others. PCA forms 
successive linear combinations of original variables that “account” for 
decreasing amounts of variation in the original data set. The problem 
of different units for different variables can be eliminated by 
Standardizing the variables with respect to the mean. We applied PCA to 
the standardized water chemistry data in an attempt to identify the 
chemical factors that could best be used to categorize surface waters 
with respect to vulnerability to acidification. 


The results of the PCA for standardized major water chemistry 
factors (Table 20) define three axes that account for 71% of the 
population variance. The variable weights are defined such that the sum 
of the weights squared equals one. Thus the weight is an index of the 
proportion of the variance accounted for by the principal component that 
is attributed to the individual variable. PC1 is dominated by the 
factors pH, alkalinity, calcium, and sulfate. These factors are related 
to acidity and buffering capacity. PC2 is dominated by chloride, 
sodium, and potassium in an inverse relationship, as these weights are 
negative. These are neutral ions that contribute to ionic strength but 
not to acidity or buffering capacity. PC3 is dominated by the metals 
aluminum and manganese. 
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Table 20. Principal Component Analysis of Standardized Maior Water 
Chem‘stry Factors of New England Lakes. Analysis Igcluded 
Only Lakes with Ion Balance of 0 + 75 yea ! 
Variable weights for 
principal components 
Variables T 2 3 
PH 0.46 0.12 -0.16 
Alkalinity 0.49 0.25 0.14 
Calcium 0.5] 0.25 0.08 
Chloride 0.16 -0.57 -0.02 
Sodium 0.17 -0.56 0.04 
Potassium 27 -0.45 0.22 
Sulfate 0.40 0.01 -0.22 
Aluminum -0.05 0.08 0.55 
Manganese 0.03 0.05 0.74 
Variance explained 34 58 71 


by each PC 








Next, the PCA scores for each lake were considered as the observed 
variables and compared to bedrock geolooy and soil class for each sample 
(Table 21). The mean values for PCl increase with increased bedrock 
buffering capacity, but PC2 and PC3 do not show consistent trends. The 
relationship of bedrock class to principal component variables was 
tested with analysis of variance. Bedrock classes 1 and 2 were not 
significantly different from each other (p = 0.10), but all other 
possible contrasts were different (p <0.01). For soil cless only PC2 
consistently increases with increasing cation exchange capacity (Table 
21). Analysis of variance shows that soi! classes SS, and SS, were 
not significantly different from each other (p = 0.614, but both were 
different from class NS (p <0.CC1). 


Based upon the groupings of chemical factors suggested by the above 
PCA, we conducted another analysis on a reduced data set. We defined 
four variables: alkalinity factor (sum of H iop, alkalinity, and Ca in 
veq 1°), salts (sum_of Cl, Na, and K in veq 1°), sulfate (S9 
concentration, »eq 1°), and metals (sum of Al and Mn, veq 1°}. We 
used oply lakes where the difference of anions and cations was 0+75 
veq 1°. The results of this analysis (Table 23) show that the first 
two PC factors account for near'y a1] (>99%) of the variance. PCl is 
dominated by alkalinity and PC2 by salts. These are the factors that 
also dominated the first two PC factors in the previous PCA. Note that 
PC1 (alkalinity) accounts for 88% of the variance. PC factors 3 and 4, 
although they contribute very little to the variance, are cominated by 
sulfate and metals respectively. These resu'ts indicate that the 
alkalinity factors -- pH, alklinity, and calcium -- are the most 
important chemice! factors in lake classification. 


Canonical Analysis 





In canonical analysis two sets of variables are transformed such 
that the correlation of the two sets is maximized. This is done by 
epplying a weighting factor to each variable in each set. By examining 
the weights applied to each variable one can ascertain which variables 
are most important in describing the variations in the data. Unlike 
principal components, however, the magnitude of the coefficients is not 
constrained, Thus, they cennot be used to directly estimate the 
proportion of canonical correlation attributable to each variable. 


Canonical analysis was used to relate physical factors such as 
bedrock and soil class, elevation, etc. to lake chemistry factors. The 
results (Table 23) show that reasonably high correlation (r, = 0.70, 

r, = 0.56) exists between physical and chemical factors. P 
cOrrelation coefficients were statistically significant for variables 1 
and 2 (Bartletts test, p <0.001). The physical factors with highest 
weights in canonical variable 1 are bedrock and soil class; the highest 
weighted chemical factors were calcium and alkalinity in the chemical 
set. Variable 2 assigned highest weights to elevation in the physical 
set and to calcium and alkalinity in the chemical set. 
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Table 21. Mean Principal Component Values for 
Bedrock and Soil Classes 





Principal component 











Bedrock class l 2 3 
] -0.89 -0.23 0.10 
2 -0.16 0.17 -0.03 
3 1.47 0.01 0.09 
4 4.68 1.56 0.52 
Soil class 
SS, 0.04 -0.88 0.09 
SS, -0.57 -0.64 -0.07 
NS 0.25 0.67 -0.01 
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Table 22. Principal nts Anulvsis of Combined Chemical Factors. 
All Units were yeq 1 “. The Factors were Defined: Alkalinity 
= Sum of H, Alkalinity, Ca; Salt = Sum of Cl, Na, K; 
Sulfate = S043 Metal = Sum of Al, Mn 





Principal component 








I 2 3 q 
Variable weight 
Alkalinity 0.999 0.016 -0.041 0.003 
Salt -0.016 0.999 -0.022 0.003 
Sulfate 0.040 0.023 0.999 -0.006 
Metal -0.002 -0.003 0.006 0.999 


% Variance 88 . 36 99.67 99.98 100.00 











Table 23. Canonical Analysis of Physical and Chemical 
Factors for New England Lakes 















Canonical variable 
1 2 3 4 















Physical factor weights 





Area -0.14 -0.03 0.14 1.03 











Hydrology type -0.22 -0.16 -0.27 0.52 
Elevation -0.38 0.63 -0.95 0.03 
Bedrock class 0.68 0.23 0.09 0.20 





Soil class 0.50 0.27 0.49 -0.20 





Disturbance code 0.35 -0.33 -0.83 0.02 





Chemical factor weights 





PH 0.22 0.26 -0.33 1.01 





Alkalinity -0.63 ~1.80 -0.96 -0.61 





Specific conductance -0.11) -0,.49 0.75 0.74 





Aluminum 0.17 0.30 -0.51 0.04 





Calcium 1.12 2.14 0.54 -0.71 





Chloride -0.3] -0.25 0.12 0.47 








Color 0.15 -0.03 0.06 -0.62 
Potassium 0.15 -0.02 -0.61 -0.41 
Magnes ium 0.24 0.17 -0.23 0.10 
Manganese -0.03 -0.06 -0.40 -0.10 
Sodium 0.34 -0.50 0.02 -0.49 
Sulfate 0.31 -0.21 0.39 -0.01 


Cannonical correlations 0.70 0.56 0.42 0.31 














Relationships between the hypothetical canonical variates and the 
variables in the set that comprise it were investigated by computing the 
correlation of each variable in the set with each canonical variate, 
following the method of Timm (1975). These results (Table 24) show 
bedrock class and disturbance code to be the most important physical 
factors in canonical variable (CV) 1. Among chemical factors, specific 
conductance, calcium, alkalinity, sulfate, magnesium, and pH were most 
important. These are the same chemical factors that dominate principal 
component 1. 


For CV2, elevation and soil type dominated physical factors and 
chloride, sodium, and potassium dominated chemical factors. The chemical 
factors are the same as those dominating principal component 2. The 
importance of bedrock and soil class in relation to lake chemistry is 
confirmed. The importance of elevation in these analyses may be the 
result of an interaction between elevation and bedrock and soil classes. 


Cluster Analysis and Discriminate Functions 

Application of cluster analysis and distriminate function analysis 
did not produce useful results. Discriminate functions could not 
classify surface waters to clusters with better than 50% accuracy. No 


rational basis for formation of clusters could be elucidated, and 
consequently these approaches were abandoned. 


MODELS OF ACIDIFICATION 
Henriksen Models 





Henriksen (1979, 1980) proposed that ionic composition of surface 
waters, corrected for marine aerosols, could be used to determine whether 
acidification had occurred. Ir unacidified waters calcium and magnesium 
are the major cations and bic bonate is the major anion. The regression 
of calcium plus magnesium on alkalinity, forced through the origin, 
produced the relationship alkalinity = 0.91 (Ca + Mg). As the ratio of 
calcium to magnesium is relatively constant over a wide range of 
concentrations, calcium can be substituted for the sum of calcium and 
magnesium. 


Henriksen (1979, 1980) hypothesized that the acidification of lakes 
was analogous to titration of a bicarbonate solution with a strong acid. 
The first stage of titration he termed bicarbonate lakes, with pH above 
5.5. The second stage was termed transition lakes, characterized by 
depletion of bicarbonate and pH fluctuations in the range of pH 5.0-5.5. 
The third stage was termed acid lakes, with no bicarbonate, excess 
hydrogen ion and ionic aluminum, and pH below 5.0. Working with surface 
water chemistry data sets from regions experiencing a gradient in sulfate 
deposition, two empirical relationships were derived that relate lake 
acidification to ionic composition. The first relates non-marine calcium 
concentration to lake pH. An empirical curve was derived based on the 
relationships of calcium and pH to alkalinity, separatinc lakes into two 
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Table 24. Correlations Between Canonical Variates and Individual 
Variables, Derived from the Standard Coefficients in 
Table 21 and the Simple Correlation Matrix 





Canonical variable 








1 2 3 
Physical factors 
Area -0.02 0.13 0.20 
Hydrology -0.18 -0.13 -0.24 
Elevation -0.19 0.90 -0.38 
Redrock class 0.81 0.27 -0.11 
Soil class 0.29 0.67 0.19 
Disturbance 0.54 -0.55 -0.55 
% Variance explained 18.3% 27.8% 9.9% 
Chemical factors 
PH 0.59 0.03 -0.23 
Alkalinity 0.74 0.08 -0.23 
Specific conductance 0.88 -0.78 -0.04 
Aluminum 0.05 0.17 -0.40 
Calcium 0.84 0.16 -0.04 
Chloride 0.09 -0.80 0.06 
Color 0.02 -0.15 0.02 
Potassium 0.35 -0.54 -0.39 
Magnesium 0.64 -0.24 -0.28 
Manganese 0.27 -0.04 -0,23 
Sodium 0.26 -0.87 -9.10 
Sulfate 0.7) 0.04 0.37 
*" Variance explained 29.6" 16.0% 5.9% 
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groups. The groups of lakes that plot below the empirical line have pHs 
in the range expected based or the calcium concentration and associated 
alkalinity. Lakes that plot above the line have pHs lower than predicted 
from the calcium - alkalinity - pH relationships for undisturbed lakes. 
This would result if addition of acid reduced the alkalinity and pH below 
that which would normally be associated with the calcium content, and 
calcium content itself was not affected. 


Application of this pH-calcium nomograg to data for New England 
lakes with non-marine calcium of 500 veq 1 * or less (Figure 38 
indicates that slightly more than half of the lakes (80 of 140; 57%) fal! 
above the empirical line and are presumably acidified to some degree. 
This is similar to the percentage of lakes in the historical pH deta set 
that are now lower in pH (see Figure 33). In our data set, pH is 
significantly related to alkalinity (Figure 39). The relationship 
differs from, the theoretical relationship. At atmospheric CO, (partial 
pressure 10 ~*~ atm), the theoretical relation is pH = 11.3 +°log 
(HCO. ). The difference in our data may be the result of higher 
partial pressure of C0. in these waters or of organic acids. 


The second empirical relationship derived by Henriksen (1982) uses 
the relationship of calcium or calcium plus magnesium concentration to 
sulfate concentration to predict lake pH. Lakes are separated into three 
groups. Bicarbonate lakes have low sulfate content, significant acid 
neutralizing capacity, and pHs generally above 5.5. Transitional lakes 
have intermediate sulfate content, have lost most acid neutralizing 
capacity, and have varying pH, but generally between 5.0 and 5.5. Acid 
lakes have no acid neutralizing capacity, high sulfate and ionic aluminum 
content, and pHs generally below 5.0. Henriksen (1980) used empirically 
derived equations equivalent to pH 4.7 and 5.3 to divide lakes into these 
three groups. 


Application of the calcium - sulfate relationships to data from New 
England lebkes with calcium plus magnesium and sulfate concentrations 
below 400 ueq 1° (Figure 40) did not agree well with measured data 
(Table 25). The model predicts pH to be much lower than was actually the 
case. One reason for this may be high background sulfate concentrations 
in New England as compared to Scandinavia. Accordingly, non-marine 
sulfate was adjusted by subtracting 50 veq 1° as the best estimate of 
background sulfate (see Figure 14) to yield net sulfate. Further, the 
lakes were divided into three pH categories pH »5.5, pH = 5.0-5.5, and 
pH <5.0, ard a linear a was run on calcium plus maonesium 
against net sulfate for lakes in the pH range 5.0-5.5. This produced the 
regression: Cg + Mg = (1.18 net sulfate) -9.9. The regression is 
significant (r° = 0.67, p = 0.0001, N= 23). This regression line was 
used to divide the lakes into two groups: pH »>5.5 and pH <5.0. The 
transitione!] lakes (pH 5.0-5.5) are plotted but are not included in the 
analysis. Application of this model to our data _—— 41) produces 
excellent results (Table 25). Acid and bicarbonate lakes are separated 
nearly completely. The transitional lakes are highly variable in pH. 
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Figure 38. Henriksen calcium-pH mode! applied to 140 New England 
surface waters having calcium concentrations 500 .eq 1 ~. 
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LOG ALKALINITY (UEQ/L) 


Figure 39, Relation of pH to log ( alkalinity in 182 New England 
lakes. The regression equation iS°pH = 4.75 + 0.87 log alkalinity, 
r” = 0.78, p = 0.0001, computed from the mean of two replicates 
and averaging surface end deep samples where taken. 
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Figure 40. Henriksen sulfate-pH mode! applied to 120 oon engigne 
lakes having calcium plus am sy concentrations 

Solid lines indicate regressions for waters of pH 4.6-4, 54.006. ) 
and 5.2-5.4 (5.3), after Henriksen (1980). 











Table 25. The Prediction of PH of New England Lakes by Henriksen's 
Nomograph and the Actual PH of these Lakes for Oricira! 
and Adjusted Sulfate Values 














~ Original Nomograph 
pl Range 
47 §.7-35.3 1B 
Number in pF range 7 14 101 
Predicted number in pH range 25 37 60 
Predicted correctly 6 5 60 
Percent predicted correctly 24% 14 100% 





Modified Nomograph 








humber in pH category 16 94 
Preaicted number in pH 16 94 
category 

Predicted correctly 14 92 
Percent predicted correctly Re° 98% 
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Figure 41. Henriksen sulfate - pH model as ip Figure 40, except 
sulfate is corrected by subtracting 5C uea 1 ° background, resulting 
in net non-marine sulfate. The equation of the regression line is 
given in the text. 
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This is reasonable considering that the lakes were sampled only once and 
lekes in this category will vary greatly in pH depending on recent 
climatological and hydrological events. 


Henriksen (1980) further proposed that the calcium - sulfate mode) 
could be used to predict pH of lakes if sulfate deposition were to 
change. This approach is based on two assumptions: 1) lake 
acidification is related to atmospheric deposition of sulfate, and 2) the 
concentrations of calcium and magnesium will not change if sulfate 
deposition changes. In southern Norway there is a steep gradient in 

sulfate deposition rate moving from south to north, and lake 
acidification is highly correlated with sulfate in precipitation (Wright 
st al. 1977; Henriksen 1979). \n New England, however, sulfate -] 

eng yt Her is relatively uniform, ranging only from 40 to 60 kg ha 

(National Atmospheric Deposition Program 1981). Attempts to 
rolees sulfate deposition to lake sulfate concentration or hydrogen ion 
concentration were not successful. However, Wright (R. Wright, Norwegian 
Institute of Water Research, personal communication) has compared 
atmospheric sulfate deposition to lake sulfate concentration for a number 
of areas in North America and found a significant correlation. These 
data sets covered a larger gradient in sulfate deposition than exists in 
New England. 


If lake calcium and magnesium content change with changes in acid 
loading, estimates of preacidification alkalinity will be in error. 
Given that ionic balance must be preserved, the addition of sulfuric acid 
to a lake may result in a decrease in bicarbonate concentration, an 
increase in cetion concentration (primarily calcium and magnesium), or 
both . Henriksen (1982) examined a number of data sets; some showed no 
change in calcium or magnesium concentration with increased acid loading 
but some showed an increase. The maximum increase in calcium plus 
magnesium concentration was equivalent to 0.4 times net non-marine 
sulfate concentration. 


Wright (personal communication) has extended the approach of 
Henriksep to a,more generalized acidificatjon equation: 
0.91 (Ca + Mg” ) - HCO, + H + Al = net SO, where the asterisk 
denotes that concentrations have been corrected for marine aerosols. The 
left side of the equation represents loss of alkalinity, or 
acidification, which is balanced by input of sulfate. The relation 
between the various components of the loss of alkalinity term will vary. 
H and Al will only be important below pH 5.0, and HCO, will only be 
important above this pH. Ca +,Mq may change, probably within the 
range of 0 to 0.4 times net SM, , but the actual change cannot be 
predicted at present. 


* * * 
For New England lakes with Ca + Mg and net SO, below 
400 veg 1“ the average calculated loss of glkalinity is to 
yeq | “, and the average calculated net SO, is 61.6 ,eq 1 
These values are not different statistically (paired to test, 


9) 











t = -1.46, * 0.15, N = 168), as expected. However, the regression of 
loss of alka inity on net sulfate for individual — pe hy; not 
significant (loss of alkalinity = 27.6 + 0.36 net SO = 0.05). 

This is pretebly the result of veriability ia lake claslatry. iedrelew, 
Climatology, etc. The data we have are not precise enough to adequately 
apply the acidification equation to an individual lake. However, the 
equation does appear to apply to the average condition in lakes in the 
region, and may be applied to predict average changes on a regional, but 
not individual lake, basis. 


To apply the acidification equation as a predictor equation the 
following relationships are ee {R. Wright, personal communication). 
The fractiog of the change in net SO, that is compensated by change 
in Ca = Yq is,given by, the sector f, which varies a 0 cad 
0.4: (Ca + Mg } = (Ca +Mq ), x F (change in net SO, ) where 
the subscripts p Bnd t denote prédicted and today, respéctively, Lakes 
above a threshold level of Ca” + Mo ) will still have bicarbonate 
and pH will be above 5.5. Lakes,below,this threshold will generally be 
below pH 5.0. The tbresho]d (Ca + Mg ) wi]l be related to net 

: threshold (Ca + Mg ) = 1.1 (met SO, ) where the, 
fattor 1.1 is derived from 1/9.91. A change in net SO,” loading will 
change the threshold Ca + Mg. 


The predicted portion of New England lakes that wil] be acidic based 
on these relationships and for various levels of net SO, are shown 
in Figure 42, In Figure 42a the factor F = 0 is used afd in Figure 42b 
the factor F = 0.4. Sulfate = 60 represents — conditions in New 
Eng lend where calculated average net SO, is 60 yeq | Sulfate , 

30 and Sulfate = 120 represent a halving and doubling of net SO, , 
respectively. For F = 0, the model predjcts that 16. oS *. = lafes will 
be acidic (pH <5) at the present net SO, of 60 ueq ; 
actual portion of lakes that were pH <5 was 8%. Thus the vote’ 
overestimates the proportion of acidic lakes. However, for F = 0.4 the 
mode! predicts that 7.5% of lakes would be acidic under present 
congitions, a value in excellent agreement with that observed. If net 
SO, is increased to an average of 120 yeq 1 “, the model predicts 
31.6% of New England lakes will 2 bate Conversely, if the average 
net SO, were to drop to 30 ,eq | the percent of lakes that are 
predicted to be acidic declines to 3.6%, 


Kramer Model 





Kramer and Tessier (1982) proposed that all cations should be 
considered in chemical acidification models, not just calcium or calcium 
plus magnesium. The theoretical relation between pH and the sum of 
cations can be calculated from equations representing carbonic acid 
weathering reactions, holding Pco, constant. Lakes that deviate from 
the theoretical line would then répresent either an input of hydrogen ion 
or an additional neutralizine substance. 
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Figure 425 Frequency distribution of lakes in intervals of 
25 weq 1 “° calcium plus magnesium. In section A a correction factor 


F = 0 is used, and in section B F = 0.4 is used. The calcium plus 
magnesium threshold between acid and bicarbonate lakes is indicated 
by the vertical line for present net sulfate concentration 

(60 veq 1°) and the projected, threshold is indicated for net sulfate 


concentrations half (30 ueq 1°“) and double (120 weq 1 “) present 


levels. Shaded bars indicate lakes that are predicted to be pH <5. 
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The application of this model to the data for New England (Figure 

43) shows that samples with high cation content and high pH plot very 
close to the line. As cation concentration decreases, the points 
increasingly deviate from the line. The deviation above the line 
indicates an excess of hydrogen ion over that predicted from carbonic 
acid weathering. Acidification was arbitrarily defined as a deviatior of 
more than 0.5 pH unit from the predicted pH. In our data, 48 of 191 
(25%) lakes and streams were classified as acidified by this 
criterion. Many of the acidic lakes fal! in this category and are 
presumably acidified as the result of atmospheric deposition. However, 
some of the most acidic lakes in our samples do not fit this model very 
well. Actual Pco, in tne sample was calculated from the equation 
pPco pH measuréd - pK + pHC0,” where p signifies the negative 

oe tnd pk 32 2.75. The measured pH was then corrected by normalizing 
re, to 10° atm., using the equation pH corregted = pH measured + 

5"pPco, where 2.5 is the negative log of 10 » and pPco 
the value calculated as described above. Some acidic sampl@s had a very 
low calculate” Pco., resulting in a very large pH correction. 
Inasmuch as Pco S not measured directly, it is not known whether 
this correction”is realistic. 


Thompson Mode! 





Thompson (1982) has proposed a cation denudation model of 
acidification, which can be expressed as a rate model or a concentration 
model. The rate mode! depicts the export of cations from a watershed in 
runoff. This model is very similar to the Henriksen model, except that 
the sum of cations replaces the sum of calcium and magnesium, and pH is 
predicted theoretically rather than empirically. The mode) allows 
prediction of pH based upon the relationship between bicarbonate ion and 
hydrogen ion in the carbonic acid system: pH = pK + pPco, - pHCO, 
In this model, pK is assumed to be 7.75 and pPco, is set at 2.5, which 
is reasonable for our samples based on the calcul$ted values of Peo, 
obtained above. The mode! assumes that the cations are balanced by 
culfate and bicarbonate. Sulfate is assumed to originate from 
atmospheric deposition (all ions are corrected for marine aerosols), and 
no correction for background sulfate is applied. 


The application of this model to our data (Figure 44) gives fair 
results. Lower pH lakes generally fit the model better thar higher pH 
lakes (Table 26). Subtraction of background sulfate does not improve the 
fit of the model significantly. This model depends on the relation of 
excess sulfate to acidification, as does the Henriksen model. As such, 
it could be used to predict pH changes for given sulfate loads if the 
relation between surface water sulfate content and atmospheric deposition 
of sulfate were known. This model assumes that cations will be mobilized 
by sulfate and therefore does not depend on stable cation concentrations. 
The cation denudation model can also be used to predict the amount of 
acidification that has already occurred. All points falling below the 
5.1 line have cations in excess of bicarbonate and are presumed to be 
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Figure 43. Application of the Kramer mode! to 19) 

New England lakes and streams. The diagonal line is the 
computed relagign for carbonic acid weathering assuming 
a Pco, of 10 ~*~ atm, 
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Figure 44, Application of the Thompson model for 142 New England lakes 


and streams. Samples with © cations greater than 400 veq 1 ~ were 
omitted. Measured mean pH is shown by the symbols. 
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Table 26. The Prediction of PH in New England Lakes and 
Streams by Thompson's Cation Denudation Model 





PH ran 
4.0-5.0 §.1-6.2 6.3-6.5 6.6-6.8 >6.8 





Measured sulfate 


Actua! number 16 60 27 20 19 
in pH range 

Number predicted 36 57 33 14 2 
in pH range 

Number predicted 14 4 6 2 l 
correctly 

Percent predicted 39% 60% 6% 14% 50% 
correctly 
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acidified. The deviation below the line is an indicator of the degree of 
acidification. On this basis, 36 of 142 (25%) of the lakes and streams 
we surveyed would be classed as acidified. This is a smaller proportion 
than that predicted by the Henriksen model, but agrees well with the 
Kramer model. 





CONCLUSIONS 


A substantial portion of the surface waters in New England 
that were not affected by direct human disturbance were found to 
be very low in acid neutralizing capacity and therefore 
vulnerable to acidification. Approximately half of the waters 
surveyed were classed as vulnerable by several! indicators, 
including alkalinity, calcium, end Calcite Saturation Index. The 
vulnerable waters were located in clusters. The locations of the 
clusters generally coincided with areas underlain by dDedrock that 
was low in acid neutralizing capacity. Both comparisons with 
historical water chemistry data and application of chemical 
models of acidification indicated that approximately 60% of the 
waters surveyed had been acidified, as indicated by a reduction 
in pH, alkalinity, or both. Presently about 8% of the waters 
surveyed were acidic (pH ¢5.0), had elevated aluminum 
concentration (>200 wg 1), or both, conditions that may be 
toxic to sensitive fish species. 
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APPENDIX A 


The data collected are presented in tabular form, subdivided 
by state. For each state the following information is presented: 


water code a unique code number assigned to 
each water sampled 

name : the name of the water 

date: date the water was sampled 

latitude and 

longitude: the location of the water in 


degrees, minute; and seconds of 
latitude and Icngitude 


elevation: elevation of water in meters 
area: the area of lakes, in hectares, 

or width, in meters, for streams 
type: 1 = lake, 2 = stream, 3 = impoundment 
hydrology: 1 = first order stream, 2 = second 


order stream, 3 = third order stream, 
4 = fourth order stream, 5 = drainage 
lake, 6 = seepace lake 


bedrock: 1 = very low in buffering capacity 
2 = moderately low, 3 = moderately high, 
4 = very high 
soil type: 1 = Sl, sensitive; 2 = S2, sensitive 
50%; 3 = SS1, slightly sensitive; 
4 = $S2, slightly sensitive 50%; 
5 = NS, mostly non-sensitive 
disturbance: 1 = more, 2 = slight, 3 = moderate, 
4 = high 
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precipitation pH: 


annual 
precipitation: 


depth: 


temperature: 


pH: 


fixed endpoint 
alkalinity: 


inflection point 
alkalinity: 


specific 
conductance: 


color: 


chloride: 


sulfate 


aluminum: 


calcium: 


mean weighted annual pH of 
precipitation as obtained from maps 
produced by the National Atmospheric 
Deposition Program 


amount, in mm, for the most recent 
year available as recorded at the 
nearest U.S. Weather Bureau station 


depth at which the sample was 
collected, in meters 


temperature of water at the time 
of collection, in degrees celcius 


pH of water sample, measured in 
replicate 


alkalinity of water sample determined 
by titration of pH 4.5, weq 1 “, 
in replicate 


alkalinity of water samp l¢ determined 
by Gran titration, weq 1 “, in 
replicate 


specific condyctance of water sample, 
uSiemens cm “ at 25°C, in replicate 


apparent color of water sample 
(unfiltered) by comparison to 
platinum-cobalt standards, color 
units, in replicate 


chloride concentration of water 
sample, mga 1“, in replicate 


sulfate concegtration of water 
sample, mg 1 “, in replicate 


aluminum conc ntration of water 


sample, vg 1 ", average of five 
measurements 


calcium concegtration of water 
sample, mg 1 ', average of five 
measurements 
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potassium: 


magnes ium: 


manganese: 


sodium: 


potassium concentration 


sample, mj 1 “, average 
measurements 


magnesium concentration 
sample, mg 1 “, average 
measurements 


manganese concentration 
sample, mg 1 “, average 
measurements 


sodium concentration of 


sample, mg 1 “, average 
measurements 


ae 


of water 
of five 


of water 
of five 


of water 
of five 


water 
of five 
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om 0. ? 9.1 4.1 5.9 “4 1.5? 0.44 o.99 9.C4#9 1. ».5 
101s 0.6 9.6 *7 4.6 " 4.47 0.98 1.w 0.70" 1.75 *.? 
1014 0.4 0.5 7.9 7.7 71 5.54 1.07 1.48 6.91% 1.4% i.* 
1017 9.6 0.4% 5.3 5.1 ™” 4.0% 0.64 0. 7” 0.091 1.14 ?.% 
191A 1.8 13.4 7? 7.0 h) 10.51 1.5? 7“ 0.067 7,77 1.* 
1019 0.4 1.9 4 %.? 377 0.93 0.35 .“ qa.nen 9.8% “.? 
1076 1.0 0.7 *.5 4.4 135 1.01! 9.56 0.18 0.07% 0.78 >.*% 
-_~ 1076 1.1 0.9 s.1 4.3 s17 1.17 0.99 0.19 9.07? 0. 40 he 
= 1989 ?.? ?.1 6.9 4.1 46 24.68 1.63 ?.89 0.194 1.9% ?.% 
1089 1.9 ?.? 5.4 $.0 s9 74.51 1.46? 7. %% 0.175 1.54 *.5 
1090 0.0 0.0 6.6 6.5 i 17.5? 0.31 9.58 0.071 0.4” 1.9 
1090 0o.% 9.4 6.5 6.3 i7 17.59 9.44 oO." 9.071 5.67 1.9 
1991 6.0 0.0 6.0 5.6 as 27.70 0.5? 0. Fl 90.0%%4 0."5 4.8 
1091 0.0 0.0 5.8 6.0 7 2.71 o.% 0.79 0.0%8 oO.nF 4.4 
109? 0.0 0.0 5 4.6 5 1.8#9 0.28 1.26 0.074 0.74 1.9 
199? 0.2 0.5 $.7 5.3 5 7.04 0.78 1.79 0.068 %. 34 1.” 
1093 0.7? 0.1 6.5 7.4 ? 1.7? 9.26 0.451 0.05” 0.51 .? 
1093 0.4 0.3 6.1 5.9 358 1.77 0.77 9.51 0.058 9.51 %.? 
19% a.4 ®.? °.4 9.4 ” m9, 7A 1.4% 14.44 0.095 $.1 1.9 
1094 9.1 9.1 %.5 9.4 ? 0.69 1.79 16.74 0.045 4.75 1.9 
1095 0.0 0.0 6.9 5.9 1° 4.62 0.45 9.41 0.0"? O.A8 4.4 
itil 0.4 0.5 7.0 5.5 55 1.55 916 0, %? 0.0% 0.44 *.5 
titt 0.5 0.6 4.8 5.1 50 1.5% 9.14 9,7? 0.0% 0.46 Pe | 
1112 0.0 0.0 14.1 14.4 i8 17.25 oOo. 78 0.79099 1.98 ».? 
iui? 9.0 .0 14.1 14.2 13 17.3% 9.55 4%. 70 0.01% 1.07 >? 
1113 0.3 0.2 9.6 9.0 0 74.41 0,86 1.™4 0.%11 0.75 1.9 
1113 0.0 0.? %1 4.8 73 74.37 0.48 1. 0.919 77 1.9 
11s 0.5 9.6 1.? 4.8 0 14.68 oo.’ ?.u7 0.088 0.75 2% 
ris 0.5 0.2 6.5 4.5 0 14.6! 9.“ r».14 0.c4l .7* he) 
1115 1.9 1.4 7.7 7.5 ™% 17.58 0.% “ 0.018 1.77 1.9 
1115 7.8 7.7 7.7 1.? 0 19.40 0.53 77 6.c07 1.*4 7 
114 7.7 ?.7 4.2? 7.4 ta] 36.9* 0.69 4,79 0.0%4 1.46 ?.% 
ite 0.8 0.7 7.8 4.3 0 39.2% 0.49 4.77 0.94) 1.47 7.5 
ai? i? 1.7 146.3 17.5 7? 9.19 0.4" 5.07 0.1% 1.75 *.5 
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APPENDIX 8 


The his »rical water chemistry data located for lakes surveyed 
in this -tudy are presented in tabular form: 


code: 


method: 


alkalinity: 


method: 


water code number from Appendix A 
name of lake 


last two digits of the year the 
historical data were obtained 


historical pH as reported 


1 = pH meter, 2 = colorimetric 
indicator 


historical alkalinity 


1 = fixed endpoint titration, 
2 = inflection point titration 
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